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Abstract: We study some gravitational instanton solutions that offer a natural realization 
of the spontaneous creation of inflationary universes in the brane world context in string 
theory. Decoherence due to couplings of higher (perturbative) modes of the metric as well 
as matter fields modifies the Hartle-Hawking wavefunction for de Sitter space. General- 
izing this new wavefunction to be used in string theory, we propose a principle in string 
theory that hopefully will lead us to the particular vacuum we live in, thus avoiding the 
anthropic principle. As an illustration of this idea, we give a phenomenological analysis of 
the probability of quantum tunneling to various stringy vacua. We find that the preferred 
tunneling is to an inflationary universe (like our early universe), not to a universe with a 
S^ very small cosmological constant (i.e., like today's universe) and not to a 10-dimensional (or 

a higher dimensional supercritical) uncompactified de Sitter universe. Some solutions are 
interesting as they offer a cosmological mechanism for the stabilization of extra dimensions 
during the inflationary epoch. 
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1. Introduction 

String theory is the only known candidate for a unified theory of fundamental physics. 
Recent understandings in string theory and its compactifications [1-5] have led to the 
realization that string theory has many vacua (here we include metastable vacua with 
lifetimes comparable to or larger than the age of our universe) . The number of such vacua, 
if not infinite, may be as huge as 10 100 or larger [6-12]. How we end up in the particular 



vacuum we are in, that is, the particular site in this vast cosmic landscape, is a very 
important but highly non-trivial question. One may simply give up on this question by 
invoking the "anthropic principle". More positively, one may take the optimistic view 
that there exists a principle which tells us why we end up in the particular string vacuum 
we are in. In this paper, we propose such a principle. A better understanding of string 
theory and gravity may eventually allow us to check the validity of the idea, or maybe to 
improve on it. In the meantime, we may use this proposal as a working hypothesis and a 
phenomenological tool. As a minimum, even if our specific proposal turns out to be not 
quite correct, we hope it convinces some readers that such a principle does exist, allowing 
us to bypass the "anthropic principle" . 

Since observational evidence [13] of an inflationary epoch [14-16] is very strong, we 
suggest that the selection of our particular vacuum state follows from the evolution of 
the inflationary epoch. That is, our particular vacuum site in the cosmic landscape must 
be at the end of a road that an inflationary universe will naturally follow. Any vacuum 
state that cannot be reached by (or connected to) an inflationary stage can be ignored 
in the search of candidate vacua. That is, the issue of the selection of our vacuum state 
becomes the question on the selection of an inflationary universe, or the selection of an 
original universe that eventually evolves to an inflationary universe, which then evolves 
to our universe today. Let us call this the "Selection of the Original Universe Principle" 
or SOUP for short. The landscape of inflationary states/universes should be much better 
under control, since the inflationary scale is rather close to the string scale. Here we propose 
that, by analyzing all known string vacua and string inflationary scenarios, one may be 
able to phenomenologically pin down SOUP, or at least discover some properties of such a 
principle, which may then help in the derivation of SOUP in string theory. The key tool 
we shall use here is a modification of the Hartle-Hawking wavefunction [17]. 

Even if we understand inflation completely, its density perturbations, and all aspects 
of astrophysics related to galaxy and star formations, no one should expect us to be able to 
calculate from first principle the masses of our sun and our earth and why our moon has its 
particular mass. On the other hand, we do understand why the mass of our sun is not much 
bigger/smaller than its measured value. That is, our planetary system is a typical system 
that is expected based on our present knowledge. As theorists, we are comfortable with 
this situation. Along this line of thinking, one should feel content if the SOUP can show 
that our universe is among a generic set of preferred vacua, even if one fails to show why 
we must inevitably end in the precise vacuum state we are in. It is along this outwardly 
less ambitious, but probably ultimately more scientifically justified, direction that we are 
trying to reach. 

Modern cosmology aims to describe how our universe has evolved to its present state 
from a certain initial state. An appealing scenario, due to Vilenkin [18, 19], Hartle and 
Hawking [17] and others [20], proposes that the inflationary universe was created by the 
quantum tunneling from "nothing" , that is, a state of no classical space-time. This quantum 
tunneling from nothing, or the ultimate free lunch, should be dictated by the laws of 
physics, and it avoids the singularity problem that would have appeared if one naively 
extrapolates the big bang epoch backwards in time. It is natural to extend this approach 
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Figure 1: Starting with nothing, quantum tunneling happens via a S D instanton with radius 
a = 1/H to a D-dimcnsional de Sitter universe, which then grows to a and beyond. 
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Figure 2: AS 4 xM instanton tunneling to the 4-dimcnsional de Sitter universe with a cosmolog- 
ically stable 6 spatial dimensional space M. Some examples are M = S 6 , S 2 x S 2 x S 2 , three- fold 
Calabi-Yau manifold with fluxes etc. 



to the brane world scenario, to see how brane inflation [21-26], natural in superstring 
theory, may emerge. Our understanding of superstring theory has advanced considerably 
in recent years so that it is meaningful to address this question. In this paper, we would like 
to consider some simple models motivated by superstring theory, to see what new features 
and issues may arise. We consider this as a first step towards the goal of understanding 
what string theory is trying to tell us about the origin of our Universe. 

For SOUP in this free lunch, we need to determine the relative probability amplitude 
of tunneling from nothing to a variety of universes. This then allows us to select the 
universe with the largest probability amplitude. Suppose that the probability amplitude 



of this tunneling to a de Sitter universe is given by the Hartle-Hawking wavefunction of 
the universe with a cosmological constant A, in terms of the Euclidean action Se as 

■®hh ~ e~ SE = e 37r/2GivA (1.1) 

Note that Se is unbounded from below [27] and small A is exponentially preferred. As- 
suming the dark energy observed today is due to a very small cosmological constant h-todayi 
then it is exponentially more likely to tunnel directly to today's universe (or our universe 
many billions years in the future) than to an inflationary universe 13.7 billions years ago. If 
there exists a solution with an even smaller cosmological constant, then that universe will 
be exponentially more preferred. Clearly this is worse than the naive anthropic principle, 
and the above formula must be modified. In fact, it was pointed out that the above prob- 
ability amplitude is unstable to corrections [28]. Intuitively, tunneling directly to today's 

— 1/2 

universe with its immense size (a ~ A to / ) mus t be suppressed. 

To tunnel from nothing to an inflationary universe that describes our early universe, we 
need a reason that selects a tunneling to some intermediate value of A to d ay < < A < < G^ 1 . 
(The observational data suggests that the tunneling to an approximate de Sitter universe 
with a A relatively close to the GUT scale is preferred.) As a phenomenological ansatz, 
we need to find an improved wavefunction by modifying the Hartle-Hawking wavefunction. 
Physically, we see that the above \P hh has not included the effects of matter fields and 
the gravitational perturbative modes around the de Sitter metric. Naively, one may think 
that those effects are small. Here we shall argue otherwise. 

We shall present 3 different (though ultimately equivalent) arguments for the necessity 
of such a modification: (i) destructive interference due to small fluctuations of large phases, 
(ii) quantum decoherence and (iii) space-like brane in string theory. The first argument 
is mostly intuitive, while the last two suggest the specific way the wavefunction should be 
improved. 

Our first argument goes as follows: (1) In tunneling to a string vacuum state which 
then evolves classically involves first Euclidean and then Lorentzian time. Treating time 
as a coordinate, this means complex metric (more precisely, complex lapse function) is 
involved. This implies that, in the sum over paths in the evaluation of the wavefunction, 
we should include paths with complex metric, not just real metric. This was suggested 
by Halliwell and Louko [29] and others for 4-dimensional gravity. (2) In summing over 
paths, the steepest descent method is employed. This is standard practice. Here, we note 
that there is a large degeneracy in paths. That is, very different paths yields the same 
action with an imaginary part, S = Sr + iSi. We expect this degeneracy to be lifted by 
the presence of gravitational and matter modes interacting with the classical metric. (3) 
When the phase Si is very large (i.e., exponentially large compared to n) and S fluctuates, 
the path degeneracy is lifted and the sum over paths will in general lead to destructive 
interference, so that quantum effects (quantum tunneling here) will be suppressed. This 
is analogous to the situation of a macroscopic particle (say a marble or a billiard ball) in 
quantum mechanics. Here we make the assertion that this decoherence takes place at the 
lifting of the above path degeneracy with very large £/s. We find that a large Si phase 
appears generically in the quantum tunneling to a vacuum state with 



(a) a large inflating volume, i.e., de Sitter size, 

(b) a small cosmological constant A, and/or 

(c) a large compactified extra dimensions. 

(4) In particular, individual paths in the 4-dimensional de Sitter case contribute to ^f an 
imaginary phase ~ 1/G/v A. So we expect destructive interference to suppress the tunneling 
probability in any of these situations. In particular, the resulting effect following from the 
sum over paths will suppress the tunneling to a small A universe. This result is very 
different from that suggested by Eq.( |l.l| ). However, a direct sum over such paths may 
be difficult. Usually, a rotation to Euclidean space allows one to sum over all paths and 
evaluate such an effect. This will lead to a new term in the wavefunction. So we are led to 
propose the following modification to the Hartle-Hawking wavefunction; 

* ~ ef, T = -S E ~ V (1.2) 

where the decoherence term T> is real positive. This wavefunction reduces to the Hartle- 
Hawking wavefunction for T> = 0. For SOUP, T> should be large in any of the three 
situations listed above; that is, the larger is the universe or the smaller is the cosmological 
constant, the larger is the value of V, so their tunneling is suppressed. 

The above argument can be made more precise in quantum decoherence. In decoher- 
ence, the classical metric (say, the cosmic scale factor a) is treated as the configuration 
variable while the perturbative modes around this metric and matter fields that couple to 
it are treated as the environment. The presence of the environment causes the quantum 
system to experience a dissipative dynamics, and the loss of quantum coherence results in 
the modification of the Hartle-Hawking wavefunction. In fact, the form of the leading term 
in T> can be easily found. A simple generalization of the known results [30-32] (which were 
used to justify the classical treatment of time) to the tunneling case in string theory yields 
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V = CV = C {^> V9 + - (L3) 



Here, c is a dimensionless constant and M s = 1/Va' is the superstring scale. Here V is the 
dimensionless "spatial volume" (measured in l s = 2tt/M s ) of the de Sitter (or any other) 
instanton, or the "area of the boundary" towards the end of tunneling. This is crudely the 
transition region from Euclidean to Lorentzian space (see Figures |l] and ||)). The origin 
of this term may be argued on physical grounds. Each mode of the environment supplies 
a mode-independent suppression factor so the resulting suppression factor is proportional 
to the total number of modes. Modes with wavelength longer than some fixed scale are 
unobservable and so should be traced over in the density matrix to yield the reduced density 
matrix. This cut-off implies that the total number of modes traced over is proportional 
to Vg, thus yielding the above T> term. A constant term in D may be absorbed into the 
normalization of \£. Quantum contributions of matter fields may also contribute to the 
prefactor V in ^ = Vexp(J-). 

Our third argument is indirect. In principle, T> should be calculable in string theory, 
that is, the T> term must have a calculable value for each of the potential string vacua, 



stable, metastable, or unstable. We may describe the tunneling as due to the presence of a 
5-brane. Since we are dealing with complex time (or Euclidean to Lorentzian time), we call 
such a space-like brane a S"-brane. A boundary term involving a S"-brane is proportional 
to Vg, leading us to the same term in T> (1.3). 



The Vg term suppresses universes with a small A and/or a large inflationary size. To 
suppress universes with very large (or uncompactified) extra dimensions, we need additional 
terms. (We do not expect T> to be simple.) Possible simple terms are Vf , V9V10 etc. In 
cases where the extra dimensions are dynamically stabilized (even only to a metastable 
state), we may use the cV term only (or a Vio term) to illustrate the approach, c may be 
a function of the vacuum state as well. It is important to determine the functional form of 
T> in string theory in a more careful analysis. For our purpose here, the above choice of T> 
is enough to explain our main points. 

In string theory, both Se and T> (and so T and *$>) should be calculable for each 
possible vacuum state. The inflationary vacuum state that ends in our today's universe 
should be the vacuum state that maximizes \& or T '. If this SOUP program is successful, 
we should be able to understand why we are where we are (why 3 large dimensions, why at 
most N = 1 supersymmetry, why SU(3) x SU(2) x U(l) etc.), thus avoiding the anthropic 
principle. 

As an illustration, we simply consider phenomenologically the above explicit term in 
T> that will provide the suppressions we are looking for and can be applied to known string 
vacua. To illustrate the SOUP idea, let us apply the above specific ansatz to the brane 
inflationary scenario proposed by KKLMMT [22]. We use Gn and the COBE density 
perturbation data to crudely fix M s , A and c, namely M s ~ 4 x 10 17 GeV, yA ~ 10 14 
GeV and c ~ 1CU 3 . Maximizing T fixes the fluxes that stabilize the Calabi-Yau manifold. 
We find that J- ~ 10 18 . Using the determined value of c, we find the value of T for the 
tunneling directly to a KKLT vacuum [2] without inflation (which has J- < 0). On the 
other hand, tunneling to a ten-dimensional de Sitter universe 5 10 , or other similar universes 
such as S 4 x S 6 , S 5 x S 5 etc. (with A suppied by D9 — L>9-brane pair) has T ~ 10 9 , largely 
independent of the value of c. In summary, the inflationary universe is much preferred 
among the vacua we examined. This is not surprising, since the functional form of T 
for a D-dimensional inflationary universe (with the remaining compactified dimensions 
stabilized) is, 

T= a - b - + (14) 

A (D-2)/2 A (D-l)/2 ^ - ^'^ 

As shown in Figure |3|, for generic constants a and b, an intermediate value of A is clearly 
preferred over a very large or a very small A. In particular, tunneling directly to a su- 
persymmetric vacuum is severely suppressed. To see which inflationary scenario or some 
other vacuum state is preferred via tunneling, a much more precise determination of ^ is 
required. Note that SOUP will not work if we use the alternative to the Hartle-Hawking 
wavefunction [33,34], where the sign of the first term becomes negative. In this case, 
large A (easily achieved with many brane-antibrane pairs) is preferred, which leads to the 
breakdown of the WKB approximation. 
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Figure 3: T as a function of A. The semi-classical approximation breaks down (indicated by the 
dashed vertical line at the right) when A is large. 



Although the details of this preliminary analysis is admittedly quite simplistic, it 
does open the possibility that a good understanding of the decoherence term may select a 
particular inflationary state which will then evolve to our today's universe after inflation. 
A more detailed analysis may allow us to determine more completely the functional form 
of T> (and the value of c). At the same time, as one examines more string theory vacua 
and inflationary solutions, one should be able to phenomenologically refine the functional 
form of T> as well. This is the basic point of SOUP. The reader may think that SOUP is 
too good to believe. We argue that the anthropic principle is too hard to believe; almost 
any scientifically motivated alternative is more desirable. 

Slow-roll inflation usually implies some form of eternal inflation [35,36]. With eternal 
inflation, it is argued that the origin of the universe issue is less pressing. Following our 
view point, this does not explain why the particular eternal inflationary universe is selected 
in the first place. If eternal inflation has happened eternally, the anthropic principle must 
be invoked in this scenario to explain why the particular inflationary universe is selected. 
On the other hand, SOUP will select the particular inflationary universe, even if it has 
the eternal inflationary property. A somewhat similar comment may be applied to the 
appearance of inflationary (or some other) universes from a "meta-universe" [37, 38] . We 
must first understand the origin of such a meta-universe. 

If successful, this program may also be used to make predictions that can be tested in 
the near future. In the brane inflationary scenario, cosmic strings are produced towards 
the end of inflation and they will leave signatures to be detected [39-46]. However, the 
type/property of the cosmic strings produced depends crucially on the particular brane 
inflationary scenario that had taken place. If SOUP can select the specific brane inflationary 
universe emerging from tunneling, it will give predictions on the type/property of the cosmic 
strings to be detected. 

It was pointed out that supercritical (D > 10) string vacua [47-49] may also provide 
a realistic description of today's universe. However, supersymmetry is typically absent 
in such supercritical vacua, so low energy supersymmetry will not be detected in LHC 
and similar experiments. It is fair to ask if SOUP can distinguish between critical and 



supercritical string vacua. Although SOUP clearly prefers S 10 over S for D > 10, with 
J-(D > 10) < 10 5 , it is not so clear for inflationary scenarios in critical versus supercritical 
string vacua, assuming viable supercritical inflationary scenarios can be realized. If so, 
the determination of the string coupling g s among various inflationary scenario via the 
maximization of J- may be able to distinguish them. Clearly a better understanding of T 
should be valuable as well. 

In this speculative paper, the following issues are discussed. In Sec. 2, we argue for 
the destructive interference in the path sum and propose the wavefunction that includes 
the decoherence term. We discuss the decoherence approach to find that the form of the 
leading term in D is the above Vg term. In Sec. 3, we apply this wavefunction for the 
quantum tunneling from nothing via a 10-dimensional de Sitter instanton S 10 to a 10- 
dimensional de Sitter space. The cosmological constant comes from the pair creation of 
branes, the simplest being D9-D9-biane pairs in Type IIB models or non-BPS Z)9-branes 
in Type IIA models. In this paper, we shall focus on Type IIB or F theory vacua. We 
also consider other cases, such as S 4 x S 6 , S 5 x S 5 , S 4 x S 2 x S 2 x S 2 S 3 x S 3 x S 4 , 
S 4 x T 1 ^ 2 x T 3 etc.. Quantum tunneling happens at a scale below the superstring scale, 
so that semi-classical gravitational analysis should be valid. Among this set of vacua, we 
find that, independent of the value of c, the symmetric S 10 de Sitter universe has the 
largest tunneling probability, with T ~ 10 9 . If the tunneling is via a S 4 x M with a pair of 
-D9-brane--D9-brane, only the 3 spatial dimensions in S 4 will grow exponentially (inflate), 
while the other dimensions are cosmologically stabilized, i.e., time- independent (though the 
amount of inflation is negligible due to the fast tachyon-rolling) . Introducing the dilaton 
does not change these qualitative features. 

Even if we assume that the moduli are all stabilized in today's universe, there is no 
guarantee that they should be stabilized during the inflationary epoch. (See Figure H[) In 
KKLT, the Kahler moduli are only metastabilized; that is, if they start out with relatively 
large values, they will continue to grow and the universe decompactifies to a higher dimen- 
sional universe. (In the simple case of a single Kahler modulus, the universe decompactifies 
to 10 dimensions.) This runaway solution is always present in superstring theory. So, if 
the extra dimensions grow during inflation, they may easily grow by a big enough factor 
and decompactification takes place. The presence of cosmological stabilization of the extra 
dimensions ensures that the Kahler moduli are not growing, so towards the end of inflation, 
they can fall down to the locally stabilized values, as indicated in Fig. ||. In our scenario, 
universes with large sizes of the extra dimensions are suppressed by the decoherence term. 
After tunneling, they remain unchanged during inflation. Their stabilization can happen 
before, during or immediately after inflation. 

We then consider the brane inflationary scenario [21] as realized in string theory by 
Kachru et. al. [22]. When we apply the wavefunction \t to this scenario, it actually selects 
a specific inflationary scenario with the fluxes determined. That is, after inflation, the 
universe settles down to a specific vacuum state, where all the fluxes (say (Mi,K\)) and 
the moduli are fixed. One may consider tunneling directly to such a KKLT vacuum state 
without going through inflationary stage. Such a preferred vacuum state will have different 
fluxes (say (Mq, Kq)). We find that tunneling to such a vacuum state without going through 
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Figure 4: The potential V(a) as a function of the size a of the compactified dimensions. If the size 
of the extra dimensions grows substantially in the early universe, the universe will end up with 10 
uncompactificd dimensions. So cosmological stabilization may be needed in addition to dynamical 
moduli stabilization. 



inflation is suppressed relative to the tunneling to an inflationary vacuum state which then 
ends up in the (M±,Ki) state. In the search of possible vacuum state for us to live in, this 
implies that the (Mq,Kq) vacuum state is not preferred. Of course, subsequent rolling [50] 
and tunneling from the {M\,K\) state to another [Mi-,K-i) state are still possible [51,52]. 
After a discussion on supercritical string vacua, we conclude with some remarks. 

Our argument for the decoherence term is admittedly crude and naive. However, fur- 
ther refinement of the wavefunction U 7 is likely to retain some of these qualitative properties. 
Examination of more possible string vacua will help us to pin down the properties of ^ . 
This approach offers the hope that SOUP will select a particular inflationary universe that 
evolves to a specific vacuum and thus the anthropic principle may be avoided. 

Many of the details are relegated to the appendices. 



2. Wavefunction of the Universe and Tunneling Probability 

The qualitative picture of quantum tunneling from nothing to an inflationary universe is 
relatively well accepted. However, the correct wavefunction of the universe, or equivalently, 
the formula for the tunneling probability, has remained controversial. Here we give a brief 
introduction (see Appendix A for a sketchy review). We then present our proposal on the 
improved wavefunction and the corresponding tunneling probability. We first discuss the 
issues for a de Sitter instanton in 4 dimensions. We then generalize the discussion to 10 
dimensions. Finally, we propose a phenomenological realization of the decoherence effect 
to be incorporated into the wavefunction. 



2.1 Birth of de Sitter Universes 

Let us first review the creation of an inflationary universe in 4 dimensions. (The gener- 
alizations to other dimensions are straightforward and will be discussed later.) We start 
with the action for a closed universe M, 



S = 



-±- [ # X y/\g\(R-2A) + -^; / d 3 x V / ]h\K+ I d^vfenWO (2-1) 
16vrG J M 8vrG J dM J M 



where A > is the cosmological constant and L m is the Lagrangian for all matter fields 
labeled as <ps. Here hij is the 3-metric on the boundary dM and K is the trace of the 
second fundamental form of the boundary. The extrinsic curvature K in the York-Gibbons- 
Hawking surface term will be important for our discussion. Here the universe is assumed to 
be closed, homogeneous and isotropic. To get the picture of the tunneling, consider Figure 
1, with metric ds 2 = N(r) 2 dT 2 + a(r) 2 dQ 2 , where the scale factor a(t) in the de Sitter 
region with lapse function N = —i (where t = Nt = —it) is the solution of the Einstein 
evolution equation : 

a 2 + 1 = H 2 a 2 , (2.2) 

where H is the Hubble parameter, H 2 = A/3. This describes the de Sitter space : 

a(t) = — cosh(Ht) (2.3) 

in the absence of matter field excitations. For positive t, this describes the inflationary 
universe. The Euclidean version of the action (2.1) can be obtained by replacing t — ► —it. 
This gives the Euclidean version of Eq. fl2.2|) : 



which gives the S 4 solution: 



h 2 + 1 = H 2 a 2 (2.4) 



a(T) = _Lcos(Fr) (2.5) 

This is the well-known de Sitter instanton, which is a compact space. It is half a four- 
sphere and is defined only for \t\ < ir/2H. This instanton is interpreted as the tunneling 
from nothing (i.e., no classical space-time) to an inflationary universe at r = with size 
a = 1/H and d = 0, as shown in Figure 1. We expect the vacuum energy to come from an 
inflationary potential V(4>), with its (local) maximum at A. 

The curvature of the S* 4 instanton solution is R = 12H 2 . For GjyA << 1, semiclassical 
approximations should be valid and we may calculate the tunneling probability following 
standard field theory procedures: Phh — e~ SE where Se is the Euclidean action of half a 
S instanton. However, 

which is negative. (The entropy S en t r0 py = —Se-) In fact, the Euclidean action for gravity 
is unbounded from below. This requires a closer look at the tunneling probability. 



10 



Following Hawking [53], let us start with the wavefunction of the universe *$> in the 
path integral formalism. For spatially closed universes, one may express Vl/f/iij] as 

^[hij] = f " Dg^ u e- S ^ (2.7) 

where h^ is the 3-metric and matter fields are ignored for the moment to simplify the 
discussion. Here S is the classical Euclidean action. The functional integral is over all 
4-geometries with a space-like boundary on which the induced metric is hij and which to 
the past of that surface there is nothing (see Figure 1). 

Consider the tunneling from nothing to half a four-sphere (r = —ir/2H to To = with 
a = 1/H) and then the universe evolves classically in the inflationary epoch (r > To to 
a > 1/H, as indicated in Figure [y). The action has the value (see Appendix B): 

3vr /_ .,, Tr2 ~2 - 3 



S(a) = S R + iS I = -—(l^i((H 2 a'-l)2) . (2.8) 

We may view the real part Sr as the Euclidean action due to tunneling while the imaginary 
part Si as the phase change due to the classical evolution of the inflationary universe. 

It is more convenient to let r be a parameter or coordinate (0 < r < 1), so the 
transition from Euclidean to Lorentzian time is facilitated by the transition of the lapse 
function from real (N = 1) to pure imaginary (N = —i) values. Once we are ready to 
entertain complex N, we should include the tunneling from nothing directly to an universe 
with size a, where Ha > 1. In this case, the path integral is 

\&(a) = / dwdz[i(z,a,w)ex.-p(—S(z,a,w)) (2-9) 

Jc 

where w symbolically labels the different paths in the complex time T-plane, and z parametrizes 
the instanton (see Appendix B for details). It turns out that S(z,a,w) = S(z,a). Us- 
ing the steepest-descent method, where the steepest-descent paths are the paths with 
Si =constant, we find again the above action (|2.8[) , where the saddle points correspond to 
Ha> 1. 

Next we consider the ten dimensional case, which corresponds closer to the situation 
in string theory. The analysis for S 10 is very similar to that of the above S 4 case. Other 
cases such as the S 4 x S & (or more generally S 1+ni x S n2 ) instanton does not differ much 
from the S 4 case; although the S 4 component becomes the de Sitter space, the size of the 
S 6 is fixed, with its value dictated by the same A (see next section and Appendix E for 
more discussions). The case of S 4 x M where the compactification of M is dynamically 
stabilized is also quite similar to the above case. The case that is interestingly new (from 
our perspective) is when the size of M is not fixed. One may consider either the case where 
the size of M also varies during inflation (e.g., Kasner-like solutions), or the case where 
the size of M is static during inflation. In either case, the result we are looking for does 
not differ much. Since the latter case is more novel, we shall focus on the S 4 x M case 
where M is cosmologically but not dynamically stabilized. To be specific, let us consider 
the case of 5' 1+ni x T l S n2 x T m where T 1 is fibered over S n2 . The Einstein equation 
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is easy to solve. See details in Appendix E. To match onto the earlier discussion, let us 
consider S 4 x T 1 ^ 2 x T 3 . We find that, after tunneling, S 4 turns into de Sitter space 
with exponentially growing a(t), while both the S 2 radius and the torus radii are constant 
(see Figure 4). That is, they are cosmologically stabilized during the inflationary epoch. 
Although the S 2 radius is fixed in terms of A, the torus radii are undetermined. If they 
start out too big, after inflation they will enter the regime shown in Figure 4 where they 
will simply grow and become decompactified. 

Following a similar analysis for the above S 4 case, we find that, for the 10-dimensional 
vacuum and Ha > 1, after summing over all degenerate paths: 

S(a) =S E + iSj = -^ K (l T i((Haf - 1)1) (2.10) 

where Gio is the ten-dimensional gravitational constant, Vq is the compactification volume 
of T l S 2 x T 3 , H 2 = 7A/30 and a is the cosmic scale factor of the de Sitter space immediately 
after tunneling. We see that the imaginary part Sj of S(a) is generically large if (i) A is 
small, (ii) a is large, and/or (hi) Vq is large. Generalization to other cases is straightforward. 
For a D-dimensional inflationary universe with n dimensional compactified volume V n , one 
has, for a > 1/H, 



a D/2 V n A r „ ,, 



3 



S(a) oc —fr-^- ( 1 =F i((Hay - 1)3 j (2.11) 

2.2 Lifting the Feynman Path Degeneracy 

Following the decoherence criterion, tunneling in each of the above 3 situations should be 
suppressed when the the effect of matter fields is included. Note that the above formula 
seems to suggest that, for Ha = 1, 5/ = even for large Vq and 1/A. It is our assertion 
that the lifting of the path degeneracy takes place in a way such that generic Si is small 
only if each of the factor is not big. This may be a very strong assumption. 

Introducing T = Nt, the integral over T in the action involves an analytic function in 
which the contour in the complex T-plane can take a variety of paths, a subset of which is 
shown in Figure |g. It is convenient to shift the r coordinate t ^ t + tv/2H, so r is always 
positive. For fixed z, let us consider 4 paths (among infinitely many) shown in Figure |B] to 
illustrate some of the points we would like to make. 

• w = tt/2H. This integration of T along the real axis (i.e., Euclidean time : r from 
to 1 with ./V = tt/2H) corresponds to tunneling to half a four-sphere, and then 
the de Sitter universe evolves classically as T goes in the imaginary direction (i.e., 
Lorentzian time) from tt/2H to (vr/2 =p iu)/H, where coshn = aH. This is the only 
path included in the Hartle-Hawking wavefunction. 

• w = 0. The integration is along the T = Nt path for 

"7T 



N 



iu 



/H 



with r going from to 1. 
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(7C/2-iu)/H 



Figure 5: The 4 paths for T = to T = (tt/2 — iu)/H discussed in the text, from left to right 
: w = 0, < w < ir/2H, w = tt/2H and w > tt/2H. Actually, there are infinitely many paths 
reachable simply by analytic continuation. Here coshw = aH. 



• < w < tt/2H. This corresponds to tunneling to less than a half four-sphere (with 
Ha < 1), and then continue tunneling to a. 

• w > tt/2H. This corresponds to tunneling to larger than a half instanton (with 
Ha < 1), and then continue tunneling back to a. 

These 4 paths (and many others related by analytic continuation in the complex T-plane) 
all have the same action given by S(z,a,w) = S(z,a) flB.9| ). The reason they have the 
same action is because of the gauge symmetry coming from diffeomorphism invariance [?] . 
Ignoring the measure and summing over paths via the steepest descent approximation 
discussed earlier and in Appendix B then yields the result ( |2.8| ) or ( 2.10 ) as given above. 
Let us make a few observations here : 



1. We are interested in an inflationary universe with growing a. Even if we tunnel to 
an inflationary universe with Ha < 1, a will rapidly grow to a such that Ha > 1. 
Summing over all paths that end at a, the real part of \P goes like — 3tt/2GA in ( |2.£j ) 
or ~ Vq/GioA in (|2.10| ), corresponding to tunneling to exactly half an instanton. This 
seems to uniquely fix the tunneling probability. This answer is same as that adopted 
by Hartle and Hawking [17] in the 4-dimensional case, though the argument is very 
different, since they do not consider paths involving complex metric. 



2. So far, the measure and the matter fields have been ignored. 



S e ff(z, a, w) = S(z, a) + S m (z, < 



ln/i(z, a,w) 



(2.12) 



Besides the the measure, which receives contributions from gauge fixing as well as 
other effects, S m includes light matter (both bosonic and fermionic) fields as well 
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as excitational modes arising from perturbing the metric. These modes/fields are 
collectively represented by (p. Intuitively, we expect the inclusion of these fields to 
lift the above path degeneracy. This will lead to an effective action S e ff(z,a,w) 
that is path-dependent. Its dependency on w will lead to important effects when the 
imaginary part Si of the action is large. A particularly interesting field to include is 
the inflaton field in a slow-roll inflationary scenario. In this case, A is actually the 
effective inflaton potential V((j)) and the de Sitter description is only approximate. 
Following the w = tt/2H path, scale-invariant density perturbation is generated 
during the classical evolution of the universe. For other paths, say the w = path, 
presumably no <f> fluctuation is produced during tunneling. In any case, whatever 
density perturbation generated during tunneling is expected to be very different from 
that generated during inflation. This illustrates the lifting of the path degeneracy. 

3. Suppose Im S(a, w) is very large; a small lifting of the path degeneracy due to 
S m (and fj,(z,a,w)) will introduce a perturbative change in S(z,a,w), which can 
still be much larger than n. This destructive interference among the paths is akin to 
decoherence, and any quantum effect will be suppressed, leaving behind only classical 
behavior. Since classically there is no tunneling, this will suppress the tunneling 
probability. Let us assume that decoherence takes place when Im S(a) is very large. 
Generically this happens in the following 3 situations : 

• Very small A. A phase variation due to S m will suppress the quantum tunneling. 
That is, in a model where A is dynamical, the quantum tunneling to a de Sitter 
universe with a very small cosmological constant is suppressed. With today's 
dark energy, (GA)" 1 ~ 10 120 . So a very small lifting of the path degeneracy that 
causes a tiny percentage change in the value on S will be sufficient to suppress 
the tunneling to a de Sitter universe with a very small cosmological constant. 

• Ha > > 1 . The lifting of path degeneracy will suppress the quantum paths that 
tunnel directly or close to a large a. This suggests that tunneling from nothing 
to a de Sitter universe with a not too large (say, close to but a little larger 
than H^ 1 ) is most likely. The universe then evolves classically (i.e., inflates) 
from such a value of a to a and beyond. This allows us to preserve the almost 
scale-invariant density perturbation generated during inflation. 

• Large Vq. Again, a phase variation due to S m will suppress its quantum tun- 
neling. That is, in a model where size of the extra dimension is dynamical, the 
quantum tunneling to a de Sitter universe with a very large compactification 
volume (or in the uncompactified limit) is suppressed. 

4. Since the real part of the action Sr suppresses tunneling to a vacuum with a very 
large A, we expect the desirable value of A to take some intermediate value below 
the Planck scale. This then justifies the semi-classical approximation we have been 
relying on. 
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2.3 An Improved Wavefunction 

In usual quantum mechanics, decoherence occurs when a macroscopic object interacts 
with its environment, e.g., a dust particle in the bath of the cosmic microwave background 
radiation in outer space. However, for the universe, there is no outside observer to play 
the role of the environment. Instead, the classical curved metric plays the role of the 
macroscopic object, while the perturbative modes of gravity as well as the matter fields 
interacting with classical gravity play the role of the environment, as explained by Kiefer 
and others [30,31]. Since this point is rather important in SOUP, we give a brief review in 
Appendix C. 

The Se term in ^ suppresses the tunneling to a large A universe. For small A, the 
universe is large, i.e., macroscopic, so we expect decoherence to suppress its tunneling (since 
there is no tunneling in classical physics). Phenomenologically we expect an inflationary 
universe with a moderate value of A (which is much larger than today's dark energy and 
smaller than the Planck scale) that evolves more or less classically before it ends with the 
hot big bang epoch. In de Sitter space or an inflationary universe, the cosmic scale size a 
plays the role of the configuration variable. Here we would like to find the leading term in 
T> due to the environment. 

A mildly path-dependent measure fi(z,a,w) and S m ((p) may allow us to realize this 
scenario via decoherence. For what values of A that this decoherence effect will begin to 
suppress the tunneling is a model-dependent question. The standard way to control a path 
integral that involves phases is to go to imaginary time. Such a calculation is beyond the 
scope of this paper. Instead, let us take a phenomenological approach here to find the 
functional form of the leading contribution to T>. Let us consider two approaches: S"-brane 
and decoherence. They yield the same functional form for the leading term in T>. 

• Although all paths in the Feynman path integral approach involve continuous metric, 
paths with discontinuous derivatives of the metric are allowed; that is, the extrinsic 
curvature K can have jumps so that its integral in the action can be finite. The 
discontinuity of K corresponds to a time-localized defect. In string theory, such 
"defects" are known as space- like branes or S-branes [54]. Since the defect here is 
located at complex time, we shall call it 5"-brane. More generally, a S"-brane has 
some structure, which is reflected in a rapidly changing (but smooth) K in some 
time region. That is, at the location of such a 5"-brane, the boundary term can be 
expressed as 

-J— ! d 9 x^/\h\K = +T S > I d 9 xyf\h\ ~ T S ,V 9 (2.13) 

167TG J dM J 

where Tgi indicates the S"-brane tension and Vg the 9-spatial volume. For the w = 
ir/2H case, Vg is simply the "boundary area" of the half-ten-sphere. One anticipates 
that the tension of such a S"-brane to be around the string scale (up to factors of 
2-7r). So we are led to the following ansatz: 

V ~expf = exp(-S E -V) (2.14) 

V ~ cV + ... = cV 9 /l 9 s + ... 
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where c is a dimensionless parameter, V is the dimensionless spatial volume, l s = 
27TVC? = 2tt/M s , and Vg is the "area of the boundary" at the transition from Eu- 
clidean to Lorentzian space (see Figures 1 and 2). That is, cV is the leading term in 
the decoherence term T>. Here, c may have a mild dependence on the matter fields, 
both open and closed string modes, as well as on the excitational modes emerging 
from perturbing the classical metric. The determination of the value of c will be a 
challenge. Sometimes it is convenient to introduce cV where V is the "spatial volume" 
of the de Sitter (or any other) instanton (as measured in string scale M 2 = 1/a'); 
that is, V = MgVg for critical string vacua. So c = (2n) 9 c. 

A constant term in T> may be absorbed into the normalization of \&. The cV term 
suppresses universes with a small A and/or a large inflationary size. To suppress 
universes with very large (or uncompactified) extra dimensions, we need higher order 
terms. Possible terms like Vf , V9V10 etc. can easily do the job, but terms like S^, 
Vg etc are not useful. Fortunately, a realistic string vacuum will have the compact- 
ification of its extra dimensions dynamically stabilized. In those realistic situations, 
we may assume the higher order terms in D to be negligible. So, with some luck, 
a phenomenological analysis keeping only the lowest order term in Eq.([L^) may be 
sufficient. 

Now we turn to the quantum decoherence approach [30,55,56], which gives the same 
leading term in T>. In quantum cosmology the total wavefunction in a Friedmann 
solution with cosmic scale a is 

ty(a,x n ) =ip (a)Y\_Xn{a,x n ) (2.15) 

n>0 

where ^ (a) is (up to a normalization factor) the Hartle-Hawking wavefunction, and 
Xn(o,,x n ), with amplitudes x n , (an orthonormal set at a given value of a) are the 
contributions from all other fields. Here the metric a is the configuration variable 
describing the semi-classical evolution of the universe while the higher multipoles of 
gravity and other matter fields interacting with gravity provides the environment. 
These fields interact with the gravitational contribution encoded in ip {a) and the 
effect of this interaction can be seen by considering the reduced density matrix, 
where the environment is traced over: 

/oo 
dx n xn(a,x n )x* n ( a 'i x n) (2.16) 

For tunneling from nothing (a' = 0) to a, we see that ^(a) ~ ps(a;a' = 0). In the 
absence of the environmental degrees of freedom, \I/ reduces to the Hartle-Hawking 
wavefunction. Let us first consider the finite a' case. In the presence of the environ- 
ment, the reduced density matrix takes the form 

Ps (a; a') ~ ^(a)^(a') [] exp (-/(a, a'){a - a') 2 ) (2.17) 

n>0 
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Contributions to f(a, a') from the tensor modes of the metric [30], massless minimally 
coupled inhomogeneous scalar fields [31], fermionic modes [32], Kaluza-Klein modes 
[57] and conformally coupled scalar fields [58] have been studied when both non-zero 
a and a' are in the classically allowed region. These works justify the validity of the 
classical time evolution of the universe. Here we use it to obtain the form of the 
leading decoherence term. 

All these modes yield very similar results. To be specific, let us focus on the massless 
scalar field in the S 4 x M case, extending the result to the under-the-barrier region. 
Each of its modes contributes (a + a') 2 /4a 2 a' 2 to f(a,a') [31], so we have 

/ /n TT / (a + a') 2 (a-a') 2 \ ( (a + a'fia- a') 2 \ , . 

p s (a; a') oc ]J exp -± ^ ) ~ exp -N { - ^ (2.18) 

n>0 ^ ' ^ ' 

where N is the total number of modes included in the environment. Taking the limit 
of infinitely many modes would reduce the Gaussian distribution to a <5-function, 
corresponding to an exact diagonalization of the density matrix. This means the 
scale factor a has been perfectly measured, which in turn implies that its momentum 
conjugate —da (or simply a) must have an infinite spread, which is highly non-classical 
(a squeezed state). For N = 0, i.e., no modes contribute, a would have an infinite 
spread, a rather unrealistic, highly non-classical situation. So, on physical grounds, 
one expects a cut-off on the modes so that the expoenent in pg (2.18) is finite. Clearly 



we need a regularization scheme. Fortunately, string theory provides a natural cut-off. 
Modes with wavelengths larger than the horizon size are not observable. Consider the 
above S 4 x M case. The higher multipoles can be found from an expansion of the full 
3-dimensional spherical harmonics (the boundary S 3 of S 4 ) and normal modes in a 
toroidal M. So n — ► (n, I, m, n\,ri2, ■■■, n&) where (n, I, m) is for S 3 and the remaining 
rii are for M. It is reasonable to only trace out modes with long wavelength A, 
(say > a/n in boundary S 3 , with a being some average of a and a' , and > L/rii, 
where Vq ~ L 6 ), we see that, N oc V^a 3 ~ Vg. Consider a field with wavenumber 
k = 27r/A = 2im/L on a torus with volume L. Since k is cut off at k max ~ 2tt/I s , 
the number of modes that contributes is A^ ~ L/l s . Generalizing this to a higher 
compactified d-dimensional volume Vd, we have A^ ~ V^/lg. For 4-dimensional de 
Sitter space in Z?-dimensional spacetime, we have 

V = cN = c(V 3 /l 3 s )(V D -4/l°- 4 ) = cVd^/I?- 1 (2.19) 

Now consider the case we are interested in, that is, the limit a' — > in 10 dimen- 
sions. In this limit, all modes are essentially frozen. The validity of a semi-classical 
description of the universe requires finite spreads in both a and a. To avoid either an 
infinitely spread in a or an infinitely spread in d, we must have iV ~ V^a 3 ~ V§aa! 2 . 
So for a' — » and a ~ H^ 1 ~ A^ 1 ' 2 , we have 

(y 6 aa' 2 )(a 2 - a' 2 ) 2 V 6 ^ 

V K ^32 " T^ - V 9 ( 2 - 2 °) 



a*a- 



A 3/2 



17 



By considering a few universes other than 5 4 x M, we see that the decoherence 
argument agrees with the £"-brane argument, that is 

tf ~ C(0) exp (S E - cV 9 /l 9 s ) (2.21) 

where ipl(0) provides the normalization. 

Note that the actual determination of T> requires knowing the full wavefunction ^f(a = 



0) Q2.15 ): this is the wavefunction when there is no classical spacetime. Since classical 



spacetime is fundamental in general relativity, \&(a = 0) may not be well-defined in 
quantum gravity. It will be interesting to see if string theory can fix \t(a = 0), or it 
has to be postulated. The parameter c actually depends on the particular vacuum in 
string theory. It is presumably a good approximation to take c to be constant for a 
subclass of vacua. 

The inclusion of D should suppress instantons with more complicated geometries, such as 
those considered in Ref [28] . Higher order terms involving higher powers of the Riemann 
tensor and field strengths are present in low-energy effective theory in string theory, so 
they should be included in Se [59,60]. In the rest of this paper, we shall ignore such terms 
and the higher order terms in T> as well as the variation of c. This is certainly sufficient 
for illustrative purposes. 

3. Toy Models : Ten Dimensional Gravitational Instantons 

To begin, we study a few instanton solutions in ten dimensions. We will be interested in 
S 10 as well as S 4 x M solutions in which four (1 time and 3 space after Wick rotation) out 
of the ten dimensions inflate while M of the remaining six spatial dimensions have static 
behavior during an early inflationary phase. These different solutions will have different 
topologies for the extra six dimensions. 

3.1 S w 

This instanton solution is a trivial generalization of the S 4 instanton. The Euclidean metric 
ansatz for S l+n is: 



dr l _ 



ds 2 = dt 2 + a{t) 2 ( -^— 2 + r 2 dn 2 n _ l} ) (3.1) 



The Euclidean Einstein equations for S +n , for any n, are: 

-(n-l)5-i(„-l)(„-2)(£-.^-A 

The bounce solution is given by a(t) = H^ 1 cos(Ht) with H 2 = 2K/n(n — 1) and R = 
n(n + l)H 2 = 2(n + l)A/(n - 1). For the case of S 10 , this becomes H 2 = A/36. By Wick 
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rotating the time axis one gets a closed 9 spatial dimensional inflating universe with the 
scale factor given by a(t) = H^ 1 cosh(Ht). 

Now, in string theory, A can come from the presence of extra branes. Starting with a 
supersymmetric vacuum, we can add N pairs of D9 — D9 branes, where the vacuum energy 
Pvacuum is just 27V times the D9-brane tension, 

M 10 
A = NA 1= 8irG wPvacuum = 16irG w N « (3.3) 

(Z7TJ g s 

where g s and M s = \^/ ja! are the string coupling constant and the string mass scale, re- 
spectively. For numerical calculations, we shall take g s = 1. The ten dimensional Newton's 
constant is given by 

8nGw = ~2Mf (3 ' 4) 

The entropy of a de Sitter space equals the Euclidean action, S entropy = —Se- In the 
presence of matter fields, matter excitation modes can be present. However, it is known 
that the pure de Sitter space saturates the entropy bound [61], so we expect the pure de 
Sitter instanton to have the largest tunneling probability. 

In string theory, the dilaton plays a crucial role in the gravity sector. A summary of 
the analysis including the dilaton can be found in Appendix D. The inclusion of the dilaton 
field does not change the qualitative features we are interested in. 

3.2 S 4 x S 6 , S n x M etc. 

We start with the 10-dimensional Euclidean metric ansatz for S 1+ni x S 712 : 

,2 _ ,2 , „/ n2 ( dr , J, n 2 \ , j,/_\2 ( d P , 2 -r>2 



ds< = dr' + a( T y {j—^ + r z dnf ni _ 1} ) + 6(r)' [j^ + p'dSl^j (3.5) 

The solution of the corresponding Einstein equation requires b to be constant (see Appendix 
E) with 

6 2 = (n 2 -l)(n 1+ n 2 -l) 

2A v ; 



ni(ni +n 2 - 1) 



For m = 3, this yields a 4-dimensional de Sitter space and a static S" 6 . That is, if we 
consider time to lie in S 4 after Wick rotation, then it describes a ten dimensional universe 
with inflation occuring in the 3 + 1 dimensions corresponding to S 4 while the remaining 6 
spatial dimensions in S e remain static. If we include a small time varying component to 
b(r), it will be rapidly inflated away. 

It is easy to generalize the above analysis to other cases. In general, for S ni x S n2 x 
S n3 x ..., when rotated to Lorentzian time, only the spatial directions in S ni which contains 
the time direction will inflate, while all the rest of the spatial directions remain static. 

We tabulate the probabilities thus calculated for the different geometries. (The list is 
incomplete and is representative only.) 



19 



Geometry 


-S E x W- g 


Volume xlO 15 


% x lO" 12 


T$ x 10 40 


s io 


3.94170 


3.92236 


1.25626 


1.75627 


S 8 xS 2 


2.88416 


2.89347 


1.27620 


1.20694 


S 7 xS 3 


3.07644 


3.11670 


1.29131 


1.23300 


S 6 xS 4 


3.17258 


3.20829 


1.31267 


1.17892 


S 5 xS 5 


3.17258 


3.26494 


1.34511 


1.07489 


S 5 x S 3 x S 2 


2.24004 


2.30866 


1.34511 


0.760062 


S 4 xS 6 


3.17258 


3.31351 


1.40018 


0.910688 


S 4 x S 3 x S 3 


2.40348 


2.49832 


1.40018 


0.686643 


S 4 xS 2 xS 2 x S 2 


1.50938 


1.59048 


1.40018 


1.07489 


S 3 x S 2 xS 2 x S 3 


1.58630 


1.73150 


1.51325 


0.335524 


S 2 x S 6 x S 2 


2.11506 


2.55074 


1.86691 


0.192098 


S 2 xS 4 x S 4 


2.26888 


2.75480 


1.86691 


0.207466 


S 2 xS 2 xS 2 xS 2 x S 2 


1.00946 


1.22435 


1.86691 


0.092207 



Table 1. The time coordinate is part of the first factor S n . The Euclidean action Se is for 
a single pair of D9 — D9 branes. The 9-dimensional volume of the boundary is measured 
in units of string scale, that is Volume= M^Vg. The value of N m is for maximizing T, the 
value of which is given in the last column. For small c = (27r) 9 c, integer N m = 1 maximizes 
J 7 , which has its value very close to — Se, not that given in the last column. 



For very small c, N m = 1 is the optimal choice (since N m is an integer); that is, one 
brane pair maximizes T '. In this case, we have J- ~ — Se ~ 10 9 , independent of the value 
of c. Let us make some comments here : 

1. In string theory, there are tunnelings to many possible expanding universes. So the 
relative tunneling probabilities become an important issue to pick out the most likely 
universes. 

2. There are at least 2 possible ways to stabilize the moduli in superstring theory: dy- 
namical versus cosmological. Cosmological stabilization has been explored in the pres- 
ence of string winding modes [62] and more recently in the presence of branes [63-66] . 
Dynamical stabilization has made substantial progress recently [2,3]. In this scenario, 
the universe is sitting at a metastable vacuum, with a tiny positive cosmological con- 
stant. As pointed out earlier, there is a true (supersymmetric) vacuum with zero 
cosmological constant where all the extra dimensions are decompactified. Even if 
we believe in dynamical stabilization, the universe must avoid an early phase where 
the extra dimensions go through a rapid growth (as illustrated in Figure 4). So cos- 
mological stabilization as presented above may be necessary even in the presence of 
eventual dynamical stabilization. 

3. The simplest instanton solution is the de Sitter S 10 instanton. Such a solution rep- 
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resents the tunneling of the universe to a space-filling D9 — D9 system in which 
all the 9 spatial dimensions inflate as the brane and the antibrane decay down the 
tachyon potential. The inflation offered by the D9 — D9 system is, however, very 
marginal, amounting only to a fraction of an e-fold. As this system decays into pairs 
of D7 — D7 and radiation, the 9 spacial dimensions undergo a radiation dominated 
expansion and this leads to a significant separation between the D7 — D7 pairs. This 
sizeable separation may lead to a significant number of e-folds in the ensuing inflation 
between the D7 — D7 system. 

4. If we use Linde/Vilenkin's wavefunction (see Appendix A for a brief review), ip ~ 
e -\ s E\ ^ e -10 ' , so it seems that (£ 2 ) 5 among the list in Table 1 is preferred in 
this case. However, an inflationary universe with large N (i.e., large cosmological 
constant) is also preferred. Of course, for large N, the semi-classical approximation 
breaks down and we lose all control of the estimate. Furthermore, the inclusion of 
the decoherence term would not help. 

Other interesting instanton solutions involve geometries of the form S m x T'S" x T k 
are discussed in Appendix E. 

3.3 An Application of the new Wavefunction 

Consider the function T '. If the cosmological constant is provided by N D9 brane-antibrane 
pairs, J- can be written as: 

Considering geometries of the form S 1+ni x S n2 \ 

T = ^7^ y i+«i V n2 ~ cMXx V n2 (3.8) 

where V± +ni and V n2 are the surface areas of S' 1+ni and S n2 , respectively. V ni is the area of 
the equator of S 1+ni . One should note that the surface areas depend on N via the Hubble 
radius, H^ 1 . To find the value of N = N m at which J- is a maximum, we just require that 
dJ-/dN = 0. This gives the following value for N: 



c 2 QA { } ni v(l + nif 



It = >r ,.T" B p.s» 



where v(n) is the "surface area" of a unit n-sphere: 

v{n) = 2vr (1+n)/2 /r((l + n)/2) (3.10) 

One can plug this value of N back in the expression for J- and get the extremized value of 

2 33 n 2 -l /2 9 Kn 2 )Kl + ni)] 9 1 
J-max - 3 18 (27r) 63^ ^ > H l [l/(ni) ]8 ^8 ^ Li ) 
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Though all geometries of type S 1+ni x S 712 have the same dependence of c (viz. 1/c 8 ), they 
have different coefficients in front of the c~ 8 term. In fact, the volume term (cM^Vg) lifts 
the degeneracy between the S m x S n and the S n x S m geometries. 

We see in Table 1 that, irrespective of the value of c, tunneling via S 10 has the largest 
probability. This is not surprising, since S 10 is the most symmetric instanton. S 10 tun- 
neling implies that the universe begins with all its 9 spatial dimensions uncompactified. 
Including the dilaton will not change this qualitative feature. Fortunately, as we shall see 
that, in string theory, the tunneling to a realistic inflationary scenario with 6 dimensions 
compactified is preferred over S 10 . 

A complete understanding of quantum gravity is lacking. However, as long as the 
quantum corrections are small, one can justify the semiclassical treatment that we have 
employed. The curvature of the instanton solutions is given by R = 5A/2 = (f-r) g s M 2 . 
The semiclassical description is reliable as long as R/M 2 ~ hNg s /Sir 2 « 1, or Ng s « 15. 
This implies that, for g s ~ 1, N ~ 1. In Table 1, we see that 

N m ~ c 2 x 10 12 (3.12) 

For the above semiclassical analysis to be meaningful, this condition requires c < 10~ 6 , or 

c=(2vr) 9 c<10 (3.13) 

which we shall assume is satisfied. For c <C 10, we have N m = 1 (since N is quantized and 
N = is ruled out). In this case, the value of T is actually very close to —Se, since the V 
term is completely negligible. As we shall see, the phenomenological analysis below puts 
c ~ 1CP 3 . In any case, the S 4 x M instanton like those in Table 1 is never preferred. 

4. Tunneling to Inflationary Universes in String Theory 

Now we turn to more realistic vacua in string theory. First, we consider tunneling to string 
vacua that initiate an inflationary phase. We then consider tunneling directly to string 
vacua that mimic our universe today (i.e., with a very small cosmological constant). To 
be specific, we shall study the model due to Giddings etc. [1] and Kachru etc. [2,22]. It 
will be important to study as many stringy vacua as possible. However, the analysis here 
should be sufficient to illustrate our approach. 

Start with a 4-fold Calabi-Yau manifold in F-theory, or, equivalently, a type IIB orien- 
tifold compactified on a 3-fold Calabi-Yau manifold with fluxes to stabilize all but the vol- 
ume modulus [1]. For large volume, supergravity provides a good description. The presence 
of Z)7-branes introduces a non-perturbative superpotential W that stabilizes the volume 
modulus in a super symmetric AdS vacuum. The introduction of l)3-branes in a warped 
type IIB background breaks supersymmetry and lifts the AdS vacuum to a metastable de 
Sitter (dS) vacuum [2] (the KKLT vacuum). To realize inflation, KKLMMT [22] introduces 
a D3--D3-brane pair, whose vacuum energy drives inflation [67-70]. The L>3-brane is fixed 
with the other Z)3-branes in the Klebanov-Strassler deformed conifold [71] and the inflaton 
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(f> is the position (relative to the D3-branes) of the D3-brane. This yields a potential for 
the mobile D3-brane in this D3-D3-hrane inflationary scenario 

V(p, <f>) = V F (K(r), W(p, <P)) + ^ + V DD (<f>) (4.1) 

where r is the physical size of the compactified volume and p the corresponding bulk mod- 
ulus. Vf(/C(7~), W(p, 4>)) is the F-term potential, where the superpotential W is expected 
to stabilize the volume modulus in an AdS supersymmetric vacuum, while the L>3-branes 
(the Dr~ 2 term) in the warped geometry breaks supersymmetry to lift the AdS vacuum 
to a de Sitter vacuum (a metastable vacuum with a very small cosmological constant and 
a lifetime larger than the age of the universe). Note that this non-perturbative interac- 
tion term leaves <fi to be essentially massless, that is, there exists a shift symmetry [23]. 
This shift symmetry is broken by the D3-D3 potential V d q(4>), which is very weak due to 
warped geometry. This inflaton potential is designed to break the shift symmetry slightly, 
so inflation can end after slow-roll. 

We shall study the wavefunction of tunneling to this type of string states. We shall 
find the maximum T for such inflationary vacuum as well as for the corresponding vacua 
without the extra pair of Z)3-Z)3-brane needed for inflation. In each case, the fluxes are 
determined by maximizing T . If we know c, this ansatz will select a specific stringy vacuum 
via quantum tunneling. Since we do not know the value of c, we shall instead use the data to 
fix it. Having fixed c, we can then evaluate T for any vacuum. Fortunately, the qualitative 
features are unchanged if we vary c by a few orders of magnitude. Before going into the 
discussions, we give a brief summary of the key results. 

• Using Gn and the COBE density perturbation data, M s and c can be fixed. We find 
c ~ 10 -3 . Maximizing T allows us to fix the values of the fluxes. 

• We find that tunneling to an KKLMMT-like inflationary universe (with T ~ 10 18 ) is 
much preferred over the S 10 universe (with J- ~ 10 9 ) discussed earlier. 

• This D3-D3 inflationary universe is also much preferred over the KKLT vacua without 
inflation (J 7 < 0). That is, direct tunneling to today's vacuum is strongly suppressed. 

• The KKLT string vacuum state that the inflationary universe ends in (after the D3- 
D3 branes have annihilated) is different from the vacuum that will be reached directly 
from tunneling. They have different fluxes. 

• More detailed analysis of the above model and investigation of other string vacua 
should help to narrow the uncertainty in c. A better understanding of \P will allow 
us to perform a more precise selection of a stringy vacuum state. 

4.1 Set-up 

The moduli stabilization follows the approach due to GKP [1] and KKLT [2], where more 
background and details can be found. Here we give an overall picture of the application to 
the wavefunction. Many details are saved for Appendix F. The complex structure moduli 
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are stabilized by the introduction of quantized fluxes wrapping around some cycles while 
the Kahler modulus is stabilized by a non-perturbative QCD-like effect. In particular, let 
M be the units of RR 3-form field strength in a 3-cycle A and —K be the units of NS- 
NS 3-form field strength in the dual 3-cycle B. Here, K and M are discrete parameters. 
Different choices of K and M correspond to different vacua. For fixed K and M, the 
introduction of an additional pair of D3 — Z)3-branes [22] raises the vacuum energy and 
inflation takes place. It is argued [23] that the slow-roll condition is generically satisfied so 
there will be enough inflation before the pair of branes collides, annihilates, (re)heats and 
starts the hot big bang epoch. 
Start with the metric 

ds 2 w = e 2u ^r Jlu/ dx^dx v + e- 2u ^g mn dy m dy n (4.2) 

where A(y) measures the initial size of the universe created via tunneling. Consider the 
simple case of a single complex modulus z. There are 2 + 2&2,i = 4 3-cycles, namely a 
pair (A,B) (dual cycles that intersect only once) and an additional pair (A',B'). The 
integrals are the periods defining the complex structure of the conifold. In particular, z is 
the complex coordinate for the cycle A 

z = f n (4.3) 

J A 

where f2 is the (3, 0)-form and that on the dual cycle is 

Q(z) = I Q, = In z + holomorphic (4-4) 

The three-cycles A and B are in the vicinity of the Klebanov-Strassler conical point where 
fluxes are turned on: 



F (3 ) = 2ttM 



1 JA 

H {3) = -2ttK. (4.5) 

B 

where Fr 3 \ is the 3-form field strength of the RR field and H< 3 \ is the 3-form field strength 
of the NS-NS field. Turning on —K' units of H^ on the B' cycle, we finally obtain the 
superpotential 

W = MQ(z) - Ktz - K'tz'(z) + Ae iap (4.6) 

where r = C(n\ + ie~^ and the last term is a non-perturbative interaction contribution due 
to the presence of D7-branes. Here, z' is a function of z which is generically non- vanishing 
at z = 0, 2/(0) ~ 1. 

Following KKLMMT we have 

A = 8.GT 3 ^(1---|) (4.7) 
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where tq and n are the positions of the stack of D3-brane and the D3-brane, respectively. 
T3 is the D3-brane tension: 

T 3 = l/(27r) 3 g s a' 2 (4.8) 

and R is the curvature radius of the AdS geometry: 

R 4 = Airg s Na' 2 . (4.9) 

Here N is the number of Z)3-branes. The D2> charge conservation requires 

N = ^F [ H {3) A F {3) = MK (4.10) 

To trust the low energy supergravity approximation, we assume N 3> 1. After compactifi- 
cation of the internal space we have 

M «=(2^)W (4U) 

where Mp l = 8ttGn and Vq is the volume of the Calabi-Yau manifold. For simplicity 
we consider Vq = a' 3 r®, where r is related to the imaginary part of the Kahler modulus 
p = b + ir 4 . 

Considering first only the imaginary part of p, we find 

-*,_»*«(£).£)«. (4.12) 

2GjvA (27r) y g s r 

Fluxes induce a large hierarchy of scale between tq and R: (%) 4 = e~ 87T9s ' 3 . We see 
that 5_g (|4.12| ) is exponentially large, due to the warp factor ^. For example, adopting 
the numbers in appendix C of KKLMMT, we find that — Se ~ 10 21 . This implies that, 
order-of-magnitude wise, T < 10 21 . 

4.2 Applying the New Wavefunction 

Now we are ready to apply our ansatz to the tunneling from nothing to the above infla- 
tionary vacuum in string theory. Since the extra dimensions are stabilized already, we 
may assume that the terms needed to suppress large compactification sizes are negligible. 
Keeping the leading term in T>, we have J- ~ e -S E -cM s v 9 ^ w \ ieve y 9 [ s th e spatial volume 
of the closed universe at the end of tunneling. 

For S 4 x CY, as long as the CY volume modulus and complex structures are fixed, 
we can use the effective four dimensional point of view such that Vg = Vq x V3, where 
Vq and V3 are the spatial volumes of the CY manifold and our three dimensional space, 
respectively. We have 

2vr .A. 3 . 

V 3 = Y (-)-*. (4.13) 



Using Eqs (fOD, ( fO| ) and ( |4.11[) with the identity V 6 = a' 3 r 6 = a' 3 p2 , we find 



/■i't' - a 2 ^ »%~-rr R \B 



V = cM»V 9 = c-^p-g s 2 (-T . (4.14) 

(2vr) 5 r 
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Maximizing ^ or T determines M to be M\ and K to be K\ . The details of the analysis 
can be found in Appendix F. For any given vacuum state, quantities like z'(0), \G(0)\, a 
and A etc. are all in principle calculable. For our purpose here, we shall simply adopt 
some generic choices of values. 

In GKP it was assumed |£/(0)| ~ 1 and z'(0) ~ 1, which we shall adopt. In this 
approximation, we find, in Appendix F, 

Mi ~ 0.4 x \A\ (4.15) 

i^ ~ 0.8 x m (l0~ 4 a 3 / 4 L4| 3 / 2 c^ 

Mpi 



(^)~ 115 xa^lil^c 1 / 2 



M s 



rsj 



3xlO- 4 a- 3 / 4 |A| 



0.1 x \A\ 

9s 

/■1-2-.1/4 

r ~ ( — ) ' 

a 
The value of T at the maximum is 

T max ~ 2 x 10- 18 a" 3 / 2 \A\ % r 2 . (4.16) 

or, expressing in terms of the COBE normalization, 8h'- 

•F^-l.exlO 4 ^ 3 . (4.17) 

For Sh ~ 2 x 10~ 5 , we find 

^ ~ 10 18 . (4.18) 

To trust the low energy supergravity approximation, we need r >> 1 and g s < 1. The 
first condition requires a << 1. For example the value chosen by KKLT, a = 0.1, results 
in r ~ 2. One may choose even smaller value of a to improve the low energy supergravity 
approximation. To satisfy the condition g s < 1, we simply choose |^4| > 10. 

We see that all the physical outputs are functions of a, \A\ and c. As an example, let 
us take a = 0.001, \A\ = 100. To get the COBE bound S H ~ 10~ 5 , we find 

c ~ 10" 10 -► c ~ 10" 3 

With above values of a and \A\, we get Mi ~ 40 and K\ ~ 12, ^ ~ 2 x 10 -3 and 
Mp( ~ 6 M s . The closeness of M s to Mp; may raise some concerns initially, but to find the 
actual physical scale of inflation one must take into account the warping effect. Including 
the effect of the warp factor, one finds 

p| ~ {^r)M s ~ 10 14 GeV (4.19) 

R 

in agreement with the commonly used energy scale of inflation. Summarizing the main 
results, we find that 
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1. Using only the Hartle-Hawking wavefunction, the flux M is fixed to a finite value 
(of order 10) while the flux K is not determined. The larger is the value of K, 
the larger is the tunneling probability. Since large K exponentially suppresses the 
vacuum energy, we end up with an universe with a very small cosmological constant, 
as expected. 

2. Using our improved wavefunction, maximizing the tunneling probability to an infla- 
tionary vacuum determines both M = M\ and K = K\. For a reasonable choice 
of c ~ 10 -3 , the preferred inflationary universe (that is, with the largest tunneling 
probability) easily fits the COBE data, with T ~ 10 18 . 

3. Using the above value for c, the 10-dimensional de Sitter universe and similar vacua 
have T ^ 10 9 . 

4. The KKLT vacua are very much like our today's universe. It is easy to estimate the 
tunneling probability to such a vacuum state, as done in Appendix F. We find that 
Fkklt < 0, so they have extremely small tunneling probabilities. Taking A to be 
today's dark energy value, we have J- ~ — 10 170 . 

5. To conclude, the universe prefers to go through an inflationary phase. This qualitative 
(order of magnitude) result is largely independent of the details (see Figure |3|). In 
this sense, we believe the picture is robust. 

6. Let K$ and Mq be the values for the vacuum with the largest tunneling proba- 
bility to go directly to an universe bypassing an inflationary phase. We find that 
(Mi,K\) y^ (Mq,Kq). This means that vacua (such as the (Mq,Kq) vacuum) with- 
out an inflationary road leading to them may be unreachable. 

7. The above analysis is very crude. Note that A and a are input parameters not de- 
termined by maximizing T ' . Also, instead of order of magnitude estimates, functions 
like z', Q etc. should be calculable. With a better wavefunction *&, we should be 
able to find the preferred inflationary universe with no free parameters. It will be 
interesting to find T for other inflationary scenarios in string theory [4, 25] . 

8. Although we find c ~ 10 -3 , the analysis is quite crude. As is clear in the analysis, 
some reasonable variation of the parameters used can easily change c by a few orders 
of magnitude. We believe c ~ 1 is quite possible. This calls for a more careful study. 



5. Supercritical versus Critical String Vacua 

Although string theory is known to have 10 (11 for M theory) as its critical dimension, 
higher-dimensional vacua are known to be possible [47-49]. In D > 10 dimensions, the low 
energy effective theory containing the massless graviton, dilaton and R-R fields is given by 
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(in Einstein frame) [47] 



Sd 



2k d < 



D-2 



d D xV ^l R _ 2 ( D - W ) e H/(D-2) 



V M </>V^ 




(5.1) 



where the sum runs over various RR fields in the theory, <j) is the dilaton and kjj stands 
for D-dimensional gravitational constant: 



2k 2 d = 16irG D = g 2 (2ir) 



D-3 M 2-D 



„2iD-2 / _ 
9s l s / 2vr 



(5.2) 



Note the presence of the cosmological constant term, due to the linear dilaton background. 
This term is absent in D = 10, allowing D = 10 flat spacetime. For D > 10, 4-dimensional 
flat, or almost flat, spacetime may be achieved in the presence of contributions from the 
many RR fluxes, branes and orientifold planes. In general there are many local dS vacua 
along with a global minimum, corresponding to <ft — ► — oo. The resulting cancellation can in 
principle yield an almost stable vacuum with a parametrically small cosmological constant 
that attains the value of the observed dark energy [48,49]. 

Such supercritical vacua in general do not have super symmetry. So the liklihood of 
low energy supersymmetry, a very important phenomenological issue, hinges crucially on 
which type of vacua nature prefers, critical or supercritical. We like to check if SOUP 
shed light on this fundamental question. A general study in SOUP program to find the 
dimensionality of the space-time is rather delicate and it may need a good knowledge of 
detailed properties of each vacuum. However, as a first test to illustrate the proposal, let 
us consider the D-dimensional instanton S D for D > 10. The main contribution for the 
cosmological constant is given above 



AOD-^V" 



(5.3) 



We would like to calculate J- for different S D for D > 10, and compare it to the value of 
J- for S 10 , the critical case. 

Using Eq,( 3.10| ) for the surface area of a unit S instanton along with Eqs.(p!2,5.3), 
we find 



V ys \2(D-W) 

I n(D-l) 1 c 

V 2(D-2)r( 



D-2 

2 



(D-l)(D-2) 

8vr 



D-l 
2 



D+l- 



(2n)*gi D " 2)/2 V 2 ( D ~ 10 ) r (#) 



(5.4) 



Using c ~ 10 3 and g s = 1 we find that J- quickly goes to zero for large D: (J 7 ~ 10 9 for 
D = 10, T si 10 4 for D = 11, T ~ 10 3 for D = 12 , T ~ 12 for D = 16 and T ~ 0.7 
for D = 20.) This is consistent with the general expectation that for large A, T — > 0, as 
is indicated in Figure 0. Adding a non-BPS brane or brane pair will simply decrease J- 
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further. As we increase D, the semi-classical approximation eventually breaks down, but 
we do not expect any change in the qualitative picture. Interestingly enough, we see that 
T for S D , D > 10 is negligible compared to the value T ~ 10 9 for S 10 , or to the KKLMMT 
inflationary model with J- ~ 10 18 . Naively, we may conclude that critical string vacua are 
preferred over supercritical string vacua. Unfortunately, such a conclusion is premature. 

It remains to be seen if some sort of brane inflationary scenario may be implemented in 
the supercritical string framework. This is a big if. Suppose a brane inflationary scenario 
similar to that discussed earlier can be implemented. Let us further assume that both 
inflationary scenarios (in critical and in supercritical strings) that agree with our early 
universe are among the vacua preferred by SOUP. Can SOUP then distinguish between 
them ? In this case, it turns out that higher correction terms to J- is required to distinguish 
them. To see why, let us start with Gn = Gd/Vd-4, where Vd-a is the volume of the 
compactified dimensions. Using Eq.( [4.13 ), a simple generalization of the earlier discussions 



gives 

3ttVd-4 _ 2ttV d _ 4 _ 3tt g| 

2G D A C 3(A/3) 3 / 2 lf~ 1 ~ 2G N A ° 487^^ (A/3) 3 / 2 [ ' ' 

For D-independent c, as suggested by the decoherence argument, there is no explicit In- 
dependence in T. With c ~ 10~ 3 , maximizing T yields vA ~ 10 14 GeV and J- ~ 10 18 , 
independent of D. To distinguish them, we may need a better knowledge of T . It is also 
important to check if a viable inflationary scenario may be incorporated into supercritical 
string vacua. In that case, maximizing J- and fitting data will determine g s and A. We 
can then compare the probabilities of tunneling to such a supercritical inflationary scenario 
and to the most favorable critical inflationary scenario. 

For example, fixing l s , c and setting dT/dA = yields T oc g~ A . With explicit 
supercritical inflationary scenarios, SOUP may be able to express its preference between 
critical and supercritical inflationary vacua. If, on the other hand, it is c instead of c that is 
L>-independent, then T oc g~ 4l (2i:)^' 2 ^ D ^ 10 ' . In this case, critical string vacua are likely to 
be preferred. It is clear that further study will be important to find the preferred vacuum. 

6. Remarks 

In this paper, we propose an improvement of the Hartle-Hawking wavefunction. The key 
difference in the new wavefunction is the inclusion of the the backreaction effects due 
to the interaction of the metric with the matter fields. These matter fields (multipole 
modes perturbing the background classical metric as well as other bosonic and fermionic 
modes) play the role of the environment and are traced over, providing a decoherence term 
that tends to suppress the tunneling to universes with small cosmological constants. As 
a consequence, intermediate values of the vacuum energy seems to be preferred. This is 
precisely what we want, since this allows the tunneling to an inflationary universe not 
unlike the one our universe has gone through. To conclude, let us make some remarks : 

• Here we like to point out that the above proposal should be further modified. Consider 
an instanton that leads to a universe with little or no inflation. That universe will 
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either recollapse, or reach only a very small size. As a result, the likelihood for us 
to be inside such a region will be extremely small compared to a universe that goes 
through an extended period of inflation. As a simple ansatz, the likelihood Ps of an 
observer ending in a particular universe S with size Vg should be 

Ps = ^W\ 2 e^ (6.1) 

where Vj is the volume of the initial universe just after tunneling, and Vs is the 
volume immediately after inflation. Here V is the prefactor that is calculable in 
principle. 

Following the Hartle-Hawking prescription, it seems that a universe with an arbi- 
trarily small cosmological constant is preferred. This corresponds to an arbitrarily 
large instanton, which corresponds to infinite (Euclidean) time for the "nucleation" 
bubble creation. This is clearly unacceptable, because it implies that the universe 
as big as today's universe is much more likely to be created directly from nothing 
than a much smaller universe, say our universe 13 billions years ago. We argue that 
there must exist a mechanism that suppresses this effect, such that a vacuum energy 
density acceptable for inflation is naturally selected. Including the coupling of matter 
fields and higher gravitational modes, we are led to a modified wavefunction, which 
includes a decoherence term. It is important to give a more complete argument for 
the wavefunction we proposed, and if additional modifications should be included. It 
will be a challenge to calculate the value of the coefficient c. 

String theory has many vacua/state with very different vacuum energies. We use 
the proposed wavefunction to compare the probability of tunneling from nothing to 
various geometries, albeit a limited set. By comparing the tunneling probabilities 
for the different geometries, we see that tunneling from nothing to an inflationary 
universe very close to what our universe went through is favored. That is, tunneling 
directly to a universe with a much smaller cosmological constant (like today's dark 
energy) is very much suppressed. Tunneling directly to a 10-dimensional de Sitter 
universe is also suppressed. 

Not all viable string vacua allow an extended period of inflation that ends with a hot 
big bang. Only those string vacua that are sitting at the end of a long inflationary 
road have any chance of being our vacuum. That is, we will not end up in a vacuum 
which is not close to a favorable inflationary road. In the search of why our particular 
vacuum is selected, this should cut down substantially the number of string vacua 
one has to take into consideration. 

So far, we have considered only a very limited set of geometries. It will be im- 
portant to check if a universe with 3 inflating dimensions is really preferred over 
all other geometries. For example, tunneling to a Randall-Sundrum brane world 
(i.e., a Horava-Witten model) [72,73] is worth a new analysis within SOUP. Clearly 
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more investigation is required with realistic string models to see if our particular 
4-dimensional universe is favored over other 4-dimensional universes. 

For SOUP to be phenomenologically meaningful, we must use the Hartle-Hawking 
wavefunction and not the alternative proposed in Ref. [33,34]. In practical terms, 
one may view this as a resolution of the 20 year old debate. 

Since it is clear that the string theory has many instanton solutions, we must enter- 
tain possibility of a collection of instantons each of which gives rise to its respective 
universe. This also suggests that, semiclassically the wavefunction should include the 
sum over the metric that includes disconnected pieces. This must allow the collision 
of different universes. In terms of string theory, presumably $ is a wavefunction in 
string field theory. Semiclassically, one may treat closed string modes, in particular 
gravity, as background fields with small fluctuations. Such fluctuations allow topol- 
ogy changes and so the interaction among the topologically distinct universes should 
be included. This is clearly interesting to investigate further. 

S entropy = —Se is the entropy of de Sitter space, maximizing the tunneling probabil- 
ity with the Hartle-Hawking wavefunction is equivalent to maximizing this entropy. 
Recall the von Neumann entropy expressed as a function of the density matrix, 

S vN = -Tr(plnp) (6.2) 

Going from pure-state to mixed-state increases S v n . It is tempting to relate the maxi- 
mizing of the tunneling probability with the modified wavefunction to the maximizing 
of the von Neumann entropy S v n for the tunneling process. 

• It is possible that such a SOUP cannot be derived from string theory as we know it, 
since the derivation of the wavefunction requires a knowledge of the full wavefunction 
\l/(a = 0,x n ) of nothing (i.e., no classical spacetime). This means that SOUP must 
be postulated. If so, it is even more important to take a phenomenological approach 
to learn about SOUP. In this case, it is unclear how to differentiate SOUP from 
the anthropic principle. One may view this as an attempt to quantify the anthropic 
principle. In this paper, we have taken a more optimistic view. Since classical 
spacetime are merely derived quantities in string theory, as opposed to the case in 
general relativity, ^(a = 0, x n ) should be a meaningful quantity in string theory, 
although a non-perturbative description may be involved. 
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A. Tunneling Amplitude 

The Feynman functional integral is over all 4-geometries with a spacelike boundary on 
which the induced metric is hij and which to the past of that surface there is nothing (see 
Figure 1). ^[/ly] satisfies the Wheeler-DeWitt equation : 

H* = (-G ijk i „ 5 x , - 3 R(h)h 1/2 + 2Ah 1/2 ) y[hij] = (A.l) 

\ ohijdhki J 

where Gijki is the metric in superspace: 

Gijki = 9 , 1/2 (h-ikhji + huhj k — hijhki) (A. 2) 

where 3 R is the scalar curvature of the intrinsic geometry of the 3-surface. For the de Sitter 
metric, TL reduces to, up to a constant factor, 

H = U-?--a + H 2 a 3 ) (A.3) 



2 V « 

where p = da. Upon quantization, p — > id/da, and the Wheeler-DeWitt equation becomes 

-d 2 



da 2 



+ U(a) *(a) = U(a) = a\l - H 2 a z ) (A.4) 



which takes the form of a 1-dimensional Schroedinger equation for a particle described 
by the coordinate a(t). This is illustrated in Figure 6. Adding matter fields \ gives the 
equation: 



ih 2 + U W + ^~r- 2e o) *(a, X) = (A.5) 

where e comes from the possibility of a matter-energy renormalization. Let 

*(o,x) = S n ^ n (a)n n (x) (A.6) 

then the above equation becomes two equations: 

1 f ~ d2 +X 2 ) u n {x) = (n + l/2)u n ( X ) (A.7) 

+ U(a)ip n ) = eV> n 



2 \d X 2 

1 f-d 2 4> n 



2 V da 2 

where e = n + 1/2 — eo- For e — ► the above equation for ipQ just becomes Eq.( [A.4| ) which 
describes the ground state of the universe. For nonzero e one obtains the excited states of 
the universe. 

Ignoring the matter fields, this reduces to Eq.( [A.4| ). The classically allowed region is 
a > 1/H, and its solutions are : 



tf ± (a) ~ e ±i K/b P( a ')^'T^/4 ^ Ai 
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Figure 6: The tunneling wavefunction. The growing and the decaying components under the 
barrier combine to form the outgoing and incoming waves on the right. If only the growing mode 
is present (the Hartle-Hawking wavefunction) , or if only the decaying mode is present (the Linde 
wavefunction), there will be both outgoing and incoming waves on the right. Vilenkin demands a 
particular combination so there is only outgoing wave. Decoherence will yield the outgoing wave as 
the classical solution. 



In the classically forbidden region, i.e., the under-barrier a < 1/H region, the solutions are 

§ ± (a)~e ± ^ 1/H|p(a ' )|c(a ' (A.9) 

The under-barrier wavefunction is a linear combination of ^+(a) and Vl/_(a), where ^+(a) 
decreases exponentially while ^-(a) grows exponentially with increasing a. There are a 
number of proposals in obtaining the tunneling probability from this wavefunction. For 
quantum tunneling, the Hartle-Hawking wavefunction is given by 



*[K 



>j\ 



[dg]e 



-S E [g] ^ e -$E 



(A.10) 



where Se is the Euclidean action of the corresponding instanton solution. This corresponds 
to the simple extension from quantum field theory. However, in gravity, Se is unbounded 
from below. In fact, for the de Sitter instanton solution, one has Eq. Ql.l| ). In terms of the 
under-barrier wavefunction, this corresponds to the dominance of the growing ^_(a). Our 
work is built on this interpretation. 

According to DeWitt [74] , the wavefunction at a = should vanish: ^ D\v( a = 0) = 0. 
DeWitt believes that a = corresponds to a cosmological singularity, so his condition 
avoids the occurrence of the singularity. This implies that ^+(a) and ^-(a) cancels each 
other at a = 0. Since they are comparable at a = 0, ^..(a) dominates around a = 1/H. 
This gives the same result as the Hartle-Hawking wavefunction. This same result is also 
obtained if one assumes a generic wavefunction around a = 0, so that the 2 wavefunctions 
\I/+(a) and ^-(a) are comparable around a = 0. 

If one assumes that only the decaying wavefunction is present at a ~ 0, then one 
obtains the following wavefunction, ^(a) — e + E , as proposed by Linde [33]. According 
to Vilenkin [34,75], ^f-(a) corresponds to an expanding universe while ^+(a) corresponds 
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to a contracting universe. Suppose one wants only an expanding universe in the classically 
allowed region a > l/H, that is, *(a > l/H) = \I/_(a). Then §+(a) and *-(o) should 
match around a = l/H, or 

*(a<l/if) = *+(a)-|*_(o) (A.11) 

In this case, the under-barrier region is dominated by VP+(a) and the tunneling wavefunction 
is ^t ~ e -37r/2GjvA_ Although this agrees with ^l for the de Sitter instanton, it is different 
in that 

y T = f[dg]e iS -» e-l 5E l (A.12) 



so for ordinary quantum field theories with bounded Euclidean action, it reduces to the 
answer we know. This alternative possibility was also mentioned in Ref. [17]. However, 
as pointed out by Rubakov [76] and others, the arrow of time is not predetermined since 
there is no outside observer. Furthermore, the introduction of scalar fields into the gravity 
theory will lead to catastrophic radiation of scalar fields (due to the sign flip of the matter 
field effective action). For our purpose, this proposal will lead to the breakdown of either 
the semiclassical approximation or to the unboundedness of the string vacuum. As pointed 
out in Ref [31], decoherence will select the growing mode as the classical evolution of the 
universe. 

The Hartle-Hawking wavefunction implies that quantum tunneling prefers a de Sitter 
universe with size a ~ 1/vA — > oo, which is clearly unphysical. This is a major reason 
in the search of an alternative wavefunction. Here we shall start with the Hartle-Hawking 
wavefunction. As discussed in the text, this wavefunction needs an important improvement 
due to the coupling/presence of other modes/fields. 

B. Feynman Paths with Complex Metric 

Here we follow closely the formalism of Halliwell and Louko [29]. Their normalization is 
slightly different from that used in the text. The metric ansatz is 

ds 2 = a 2 [Nirfdr 2 + b{r) 2 d^l) , (B.l) 

where a 2 = 2G/3ir. We choose r to be real and a (piece-wise) constant lapse function: 
N = 0. 

The Einstein-Hilbert contribution to the action (j2.1|) is 



where A = a 2 A/3 = a 2 H 2 and the last term above comes from partial integration of a 
second derivative term in the Ricci scalar, containing b. The extrinsic curvature is 

«--£>■ ™ 
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so the boundary term contribution to the action (2.1) is 



and 

1 f 1 bh 2 

S = -jf drN(-— + Xb'-b). (B.5) 

Taking the initial value of 6(0) = and the final value 6(1) = a, the solution of the action 

is 

6(r) = rsin(— ) . (B.6) 

r 

where r = 1/vA. However, we shall leave r unfixed, since we are interested in summing 
over four-spheres of arbitrary radius r. Instead of r, one may trade it for z, where 

* = 1 + -^ = 1 + cob( 7 ) (B.7) 

and 

r 2 = -r^- (B.8) 

2(2 — z) 



Inserting this ansatz in the action (|B.5|) and performing the the r integral we find 



a 2 ( Xa 2 - 4 Xa 2 
-ll-z+ + - 

D V Z Z 



S(z,a,T) = S(z,a) = ^ r [l-z+ — + ^- ) , (B.i 



and the path integral is just a single ordinary integral over z: 

\&(a) = / dwdzfi(z,a,w)ex.p(—S(z,a,w)) (B.10) 

Jc 

Note that we have introduced w to label the paths that have the same degenerate action 



( B.9 ), a fact that will be important later. 

We shall adopt the steepest-descent method, where the steepest-descent paths are the 
paths with Im(S) =constant. For Xa 2 < 1, the saddle points are 

z = l±(l-Aof)3 , (B.ll) 

which correspond to r = A -1 ' 2 and represent, respectively, less than or more than half a 
four-sphere with three-boundary with radius a. (We ignore another saddle point at z = —2, 
which apparently corresponds to negative A.) The action at these two saddle points has 
the values 

S(a 1 ) = -^-(l T (l-Xa 2 1 )l) . (B.12) 
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respectively. One may envision that N stays real as the tunneling from 6(0) = to 
6(1) = oi, at which point the universe starts evolving classically in an inflationary epoch, 
where a grow from a\ to a. That is, at r = 1, we switch from the Euclidean metric (for 
5 4 ) to a Lorentzian metric (for de Sitter space), i.e., the lapse function N goes from a 
real constant to a pure imaginary constant. To describe the tunneling of an inflationary 
universe from nothing, we must deal with values of N in the complex N plane. 

Once we are ready to entertain complex N, let us consider the the tunneling from 
nothing directly to an universe with size a, where Xa 2 > 1. In this case, the path integral 
has the saddle points at (again ignoring the saddle point z = —2) 



1 ± i(Xa 2 - 1) 



(B.13) 



at which the action has the value 



S(a) = --(l T i(\a 2 - 1)* 



(B.14) 



respectively. Here, the saddle points still correspond to r = X 1 ' 2 , but r < a. We still have 
6(0) = 0, 6(1) = a and b(0)/N = 1, but 



N = r 



TT 



(B.15) 
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a(r) = r sin ( — J cosh(ur) =F i cos ( — J sinh(«r) 



cosh(u) = A 1/2 



a 



These new solutions arise only when one performs the analysis with complex metric 

Now, let us come to the issue of degeneracy. Introducing T = Nt, the action ( |B. 
becomes 



S 



1 



N 



dT 



-»(£) + - 3 " 6 



(B.16) 



this is a complex integral over T of an analytic function in which the contour in the complex 
T plane can take a variety of paths, a subset of which is shown in Figure 5. 
It is instructive to calculate the extrinsic curvature. Using Eqs 



T . 3a/A / ir 

K = cotan r(— =F iu) 



and (BJ3) we find 
(B.17) 



For r in the interval < r < 1, K is in general a complex number. At r = 1 we find 



K 



3/ 



aa 



(Xa 2 



1)5. 



(B.18) 



C. Decoherence 



Decoherence introduced in Eq. ( |1.2[) is brought about by lifting the degeneracy of complex 
metric paths by introducing matter degrees of freedom. However, the concept of decoher- 
ence is relevant and applicable in a broader framework. Decoherence can be understood 
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in very general terms in quantum mechanics as arising due to the interaction between the 
system and the environment. The effect of decoherence has been explored in the literature 
in various quantum mechanical situations, including quantum cosmology. In any quan- 
tum mechanical problem there is a system, consisting of the part that we are interested in 
studying, and the environment, that consists of the part that is not directly relevant. The 
environment can interact with the system in such a way that leads to a loss of quantum 
mechanical features (i.e., a loss of coherence) and this process is called decoherence. This 
may be compared to the quantum Zeno effect where a continuous measurement of the 
system by the environment leads to the suppression of its tunneling [77, 78] . 

A simple example consists of an electron beam (the system) in the double slit ex- 
periment in quantum mechanics with a gas of molecules between the slits and the screen 
(the environment). The random interactions between the gas molecules and the electrons 
can lead to a loss of coherence and the interference pattern can get washed out leading 
to a decoherent pattern on the screen. The idea of decoherence has also been applied to 
quantum cosmology to study the semi-classical evolution of the universe. Although one 
would expect the universe to evolve as some superposition of various possible histories (for 
example, a quantum mechanical superposition of the growing and the decaying mode in 
the case of the de Sitter solution) , decoherence allows the universe to evolve classically. 

The decoherence term in Eq. (|l.2| ) ought to arise as a result of such interaction between 
the system (the gravitational instanton) and the environment (consisting of all other fields) 
[30-32]. Decoherence should forbid the quantum tunneling of a large (macroscopic) sized 
universe. Whereas in previous works [30, 31] the idea of decoherence has been applied to 



the outside barrier region of the potential in Eq.( A. 4) to study the classical evolution of the 
universe, we attribute the origin of the decoherence term in Eq.( |1.2[ ) to the decoherence 
mechanism in the under barrier region. So it is an extension of the previous ideas to the 
under barrier region. 

It is convenient to use the language of density matrix to study decoherence. It is also 
possible to cast the arguments using path integral language but we shall use the density 
matrix language here. Consider the total wavefunction |^ > of the form : 

n 

The corresponding pure-state density matrix is: 

P=|* >< *| =J2 C " C *rn\Sn > \K >< E m \ < S m \ (C.2) 

van 

where \S n > are states of the system and \E n > are states of the environment. In this 
state, the system and the environment are correlated with each other. We are, however, 
only interested in the state of the system, not the environment. So the environmental 
degrees of freedom can be traced over in any calculation of interest. The object of interest 
is then the reduced density matrix given by: 

PS = tr E \® >< *| = J2c n c* m < E m \E n > \S n >< S m \ (C.3) 
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If the dot products < E m \E n >= 5 mn , the off-diagonal terms of the density matrix vanish 
so there are no interference terms, and the system displays a pure classical behavior. The 
reduced density matrix then becomes diagonal: 



PS 



^M 2 |,S n ><<S„| (C.4) 



In general, < E m \E n > is only suppressed when m ^ n. Decoherence is the suppression 
of the off-diagonal terms, so the system behaves almost classically. In quantum cosmology 
the total wavefunction is 

*(a, x n ) = ip (a) Y[ Xn(a, x n ) (C.5) 

n>0 

where Y> (a) is the Hartle-Hawking wavefunction, and Xn{a,x n ) are the contributions from 
all other fields with amplitudes x n . These fields interact with the gravitational contribution 
encoded in V'o(o) and the effect of this interaction can be seen by considering the reduced 
density matrix: 

/oo 
dx n Xn(a,x n )x* n (a\ x n) 

Here p(a;a ) should be compared with the density matrix element p mn . In the absence 
of the environmental degrees of freedom, the density matrix element corresponding to the 
tunneling of the universe as derived from the Hartle-Hawking wavefunction is given by 
p(a = i? _1 ;0) = e~ E . However, tracing out the environment is expected to give the 
decoherence term leading to the wavefunction proposed in Section 2, 

tf(o) =p(a = H- 1 ;0) ~ e" SB-2 ' (C.6) 

The exact form of T> depends on the particular model that we consider. A calculation 
can be performed along the lines in Refs [30,31] with the environment given by the higher 
multipoles of the matter and geometry as obtained by Halliwell and Hawking [79]. That 
tracing out these multipoles leads to decoherence in the Lorentzian evolution of the universe 
has been shown in [30]. We are proposing that a similar procedure can be applied to the 
under barrier Euclidean/complex metric regime and that can lead to the appearance of 
the decoherence term D. In this sense, the inclusion of the complex metric paths and 
the S -brane is related to the appearance of decoherence by tracing out the environment. 
The decoherence should appear in string theory by tracing over the closed and open string 
degrees of freedom. Hopefully, string theory provides a natural cut-off that is absent in the 
analysis so far. 

Although the exact form of T> is expected to be obtainable only from a complete 
analysis in string theory, we shall give a heuristic derivation here which shall justify the 
form of T> we have claimed. Consider a geometry of the type S 4 x M where M is a 
6-dimensional compactified space of volume Vq. Consider only the tensor modes as the 
environmental degrees of freedom as in [30]. To trace over the environment, we need to 
know the wavefunction of these tensor modes in the under-the-barrier Euclidean region. 
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However, as shown in Ref. [79], the form of the tensor mode wavefunctions is the same both 
for the Euclidean and the Lorentzian regions and this is especially true for large n (small 
wavelengths) modes. In particular, for both the Euclidean and the Lorentzian cases, the 
tensor modes are given by: 

raa 2 \ 1/4 ( n .dS 2 1 2a 2 



Xn(a, x n ) = ( — I exp I -2i—x n - -ne a x n ) (C.7) 

where S is the Hamilton- Jacob i parameter. The phase parts differ for the Euclidean and the 
Lorentzian regions. But it is the amplitude part of the tensor modes that are important for 

determining the decoherence term and they are the same, viz. I ^- J exp (-jjje 20 ^) 
for both regions. Tracing over the tensor modes involves finding 

/ 2 '2\ _1 / 2 

dx n x* n (a,x n )xn(a ,x n ) = \/2 ; — (C.8) 



This leads to the result we have used in Eq.(2.18) 



p s (a;a)oc_[J_exp -^— ~ exp (-N —^ (C.9) 



4aV 2 J l V 4aV 2 

n>0 7 v 

Tracing over only long wavelength modes, the number of such modes is proportional 
to the spatial volume. So it is only natural that N should be proportional to the 9-volume 
Vg. Writing Vg as Vqo, 3 , we have argued in Sec. 2 that o 3 = aa' 2 for consistent decoherence 
to take place. Any other choice would lead to either zero or no decoherence Here we are 
staying in the Euclidean region and starting with a finite a' and a; we then take the limit 
where a' — ► 0. Decoherence is seen if we choose the cut-off scheme N ex V^aa' 2 . Strictly 
speaking, the whole calculation in quantum gravity is ill-defined at a' = where classical 
space-time does not exist. Another way to see the difficulty is that, close to a' = 0, the 
separation between lower and higher multipoles becomes less and less accurate. The entire 
perturbative calculation done in Ref. [79] breaks down at a' = 0. Since classical spacetime 
are merely derived (i.e., not fundamental) quantities in string theory, the a' = limit 
should be meaningful and hopefully calculable. There have been discussions about the 
calculation of the wavefunction of the universe in a string theory context in [80] . 

Perhaps we can convince the reader that the decoherence term Pcx Vg by considering 
the case where the universe starts off at some finite a' with energy less than the height of 
the barrier as described by a nonvanishing e parameter in Eq.(A.7). After tunnelling, it 



becomes a deSitter state on the other side of the barrier in Figure 6. The calculation of 
the decoherence term can now be done with tensor modes and the final result should have 
Doc Vg. We expect the qualitative result V oc Vg to be true for this analysis. One can then 
take the limit e — > and obtain the result for the tunneling from nothing. A macroscopic 
universe would not tunnel because of decoherence. 
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Figure 7: S 4 instanton 



D. S w with Dilaton 

In this section we investigate the impact of dilaton on the action for S 10 . The action with 
dilaton included in string frame is 



S 



1 



16vrG 



-20 



d w xyJ\g\e~ Z(p [R - 2A e* + 4(V<£)' 



(D.l) 



After the Weyl transformation g mn — > e^g mn , the equations in Einstein frame are 

,2 1 
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We are interested in the instanton solution with the following properties: 

a(h) = a(t 2 ) = 

a(t) < a(t ) 

o|i = 0|t o = ° 



(D.2) 
(D.3) 
(D.4) 

(D.5) 



The first two conditions correspond to having a closed instanton with spherical topology. 
The third condition comes from the analyticity and real-valued properties of the solution 
in transition from Euclidean space-time to the Lorentzian space-time (see Figure pi). We 
show that, with the dilaton included, there are apparent singularities at t = t\ and t = t 2 . 
In other words, £2(^1) corresponds to big-bang (big-crunch). 

To begin the proof, suppose the solution is non-singular. The point of special care are 
t\ and t 2 , where a(t) = 0. We focus our investigations near t = t 2 - The same discussion 
applies equally to the t\ case. To have a regular solution, all scalars like R, R mn R mn , 
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RmnpqR mnpq \ ■■■ constructed from the metric components and their derivatives must be 
regular. In particular, we have 



-timripq-ti — Ot) 



rr'+U'^-i 



a \ a 



(D.6) 



Clearly, R mnpq R mnpq is regular if both (|) and (|) 2 — \ are regular. Since a(t) vanishes 
at t 2 by construction, ^ will diverge such that (|) 2 — K is finite. Suppose 



a(t)~a[(t-t 2 y + j3(t-t 2 ) q ] , 0< P <q. (D.7) 



We find 



1 „ plt-tz)- 1 + Piq-pW-h)*-*- 1 

a 

- a ~p(p- \){t - t 2 )~ 2 + 15 [q(q - 1) - p(p - 1)] (t - t 2 y-?~ 2 (D.8) 

To have (|) 2 — \ and | regular, we need a = —1 and q — p> 2, respectively. This implies 
that q > 3. So 

o(t) ~ (*2 " *) + P(t - t 2 ) q , q>3. (D.9) 

On the other hand 



d, 9 1 a 
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Near <2> by using Eq(D.9), we find 

<j> 2 ~ 8/3q(3 - q)(t - t 2 ) q - 3 . (D.ll) 

Previously we concluded that g > 3, which indicates that </> — > with some positive power. 
Now consider the scalar field equation. It is more instructive to write Eq( p.4|) in the 
following form 

^-{a 9 (j)) = 3Aa 9 e¥ (D.12) 

Define F(t) = a 9 (ft. At t = 0, we demand (ft = 0, so -F(O) = 0. Near t = t 2 , the regularity 
assumption of the solution requires that a(t) ~ (t 2 — t), (ft — > 0, which implies that F{t 2 ) = 
as well. With the conditions F(0) = F(t 2 ) = 0, the mean value theorem of analysis implies 
that F vanishes at some point in the interval [0, £2]- But F = 3Aa 9 e~ is always positive 
definite. We conclude that regularity at t = t 2 is a false assumption and the system develops 
a singularity at t = t 2 . It is easy to see that the only property of the dilaton exponential 
potential which was used to verify the existence of the singularity is the fact that ^C > 0. 
This means that the singularity is a generic properties of any scalar field with monotonic 
potential. 



41 



In order to evaluate the action we need to find the form of the singularity. Using 
Eq.flD7?D, from Eq. flDlOD we find 



r~^ 



1 M J ^ w - 1 s irWr \ 



To have cp 2 > 0, one requires p < 1, which implies that <fi ~ j—^- up to a coefficient. The 
possibility of p = 1 is excluded. The previous discussion indicates that F(t) = a 9 cj) should 
not vanish at t = t 2 , which is not possible for a(t) ~ (t — t 2 ) and <ft ~ j^tT' This will fix 

the coefficient of 4> and we find <ft = ^ a-t ) anc ^ ^ ~ ^-^\fv l n l^ ~~ ^2 1- Using the scalar 
field equation, one can easily see that the positive branch solution for (p is not allowed and 
matching the most singular power of , t _ t , will fix p = g . In summary, we find 

a(t)~(f 2 -t)5 , <p{t) ~ -ln\t - t 2 \ . (D.14) 

o 



:>„-< 



To evaluate the action, from Eqs( p.2j ) and (|D.3|) , (j) and Ae 2 can be expressed in terms of 
a(t) and its derivatives. We find 



; ' d w xa 9 
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(D.15) 



Near t = t 2 , the singular terms in the bracket cancel each other and the integrand behaves 
like (t 2 — t)& and the action is finite. 

Although the analysis are carried out only for S 10 , one can expect the same qualitative 
result apply for other instantons considered in section 3 and the dilaton inclusion will not 
change the result. 

E. S 4 x S 6 and Other Instantons 

We start with the 10-dimensional Euclidean metric ansatz for S 1+ni x S n2 (ni+ri2 + l = 10) 
with the metric (|3.5|). The Einstein equation G^ v = 8ttGiqT^ c = —(Aag^ n , A^g^ 2 ) 
becomes (the time-time component, the space-space component in S 1+ni and the space- 
space component in the S U2 respectively) 

ni(m-l) (a 2 1\ n 2 (n 2 -l)(b 2 l\ ab 

^-^J 2 [tp-p)- nin2 a6^ Aa (K1) 
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a b ni(n\ — 1) (a 2 1 
-^a-^-^b 2 {*-* 

(n 2 -l)(n 2 -2) (b 2 l\ ab 
2 [¥-¥ - ni{n2 ~ 1] ab =Ab 
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where the dot is with respect to Euclidean time r. (Different As may be realized by having 
some D9 — D9 branes plus, for example, D(l + ni) — D(l + m) pairs of branes behaving like 
dust in the extra n 2 dimensions.) The instanton solution corresponds to a S 1+ni bounce 
with a static S n2 (constant 6(r) = b): 

- d 2 + 1 = H 2 a 2 (E.2) 

a( T ) = # -1 cos(#t) 
H 2 = 2(ni - l)A a - n 2 (rn - l)((A a - A b ) 
ni(ni - l)(ni + n 2 - 1) 
,2 _ (n 2 - l)(ni +n 2 - 1) 
~ (A a + A 6 )+m(A a -A 6 ) 

Here S 1 " -2 is cosmologically stabilized, with radius b. For A a = A& = A, the above result 
reduces to 



(n 2 - l)(ni +n 2 - 1) 
2A 

h»- 2A 



(E.3) 



ni(ni +n 2 - 1) 



For the case n 2 = 0, the last component of the above Einstein equation is absent and we 
have H 2 = 2A/n(n - 1). For S 10 , this becomes H 2 = A/36. 

For the case of S x S 6 with equal A, n\ = 3 and n 2 = 6, and we get b = y^O/A anc j 
i7 2 = A/12. By Wick rotating the time axis to Minkowski metric: r — ► it, we have 

a(t) = H^coshiHt) (E.4) 

this is just the usual 1 + 3 dimensional deSitter space-time with a static S e making up the 
extra dimensions. It describes a universe that is contracting at t < 0, reaches its minimum 
size (a m i n i mum = 1/H) at t = and then expands for t > 0. For large t, the expansion is 
exponentially, corresponding to an inflationary universe. Since the cosmological constant 
is due to brane-anti-brane pairs, this universe is always classically unstable. The presence 
of tachyonic modes will lead to their eventual annihilations, which will end inflation. On 
the other hand, the above Euclidean solution describes a compact space; it is defined only 
for \t\ < ir/2H. Together, they describe the quantum tunneling process (See Figures 1 and 
2). Starting from "nothing" (i.e., no space and no time), the instanton solution bounces 
at the classical turning point a = 1/H (where a = 0), and continues at Minkowski time t. 
The Euclidean action Se an d the spatial volume are easy to determine. As an example, 
consider the case of S 4 x S & . The geometrical factor J ^/gd 10 x factorizes to A^ x Aq. The 
radius of S* 4 is just the Hubble radius, H~ l . The radius of S* 6 is b = y^O/A, so 
A 4 = |vr 2 i7- 4 and A 6 = ifvr 3 & 6 . So 



h d "'*=mk- x T xA * xA > (E ' 5) 
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The generalization to S ni+1 x S n2 x ... x S nk instantons is completely straightforward. 
One finds: 

ai (t) = H^cosiHt) 



2A 

where Z = 2, 3, ..../c and D = ^ i=1 Jij + 1- The cosine dependence of S ni+1 becomes a cosh 
dependence in Lorentzian space-time and the solution describes, after tunneling, inflating 
ri\ + 1 dimensional space-time with X)j= 2 ni extra dimensions that are static and have the 
geometry of S n2 x ... x S nk . 

One may also consider instantons with torus fibered over sphere : S ni+1 x T l S n2 x T m . 
For example, consider T 1 fibered over S 712 as implied by the notation ... x T 1 S ri2 x .... The 
metric is given by: 

ds 2 = dr 2 + a(r) 2 ( 1 ^+r 2 ^ i _ i) ) (K7) 

+b 2 (^$ + P 2 ^ 2 n2 _ 1} ) + c{0fda 2 + d 2 YJIU dP? 

where a and Pi are the scale factors of the torii. 9 is the azimuthal coordinate of S 712 and 
the 6 dependence of the torus scale factor c(9) expresses the fibration of the torus T 1 on 
S n2 . c{6) and d are taken to be time independent as we are looking for solutions where the 
extra dimensions are time independent. The Euclidean Einstein equations for this metric 
are: 

ii , , / a 2 1 \ 

— + m-1 --- =0 
a \a z a z J 

-n 1 ---n 1 (n 1 -l)^-^j--(n 2 -l)(n 2 -2) p - p ^cot^j=A 

-nx- - -n l( n, - 1) ^- - -j - -(n 2 - l)(n 2 - 2)^ - - 2 (+-li j = A 

a 1 n/« 2 1 \ 1 . 1 

-m- - -Mm " i) ^ " ^J " ^ " % = A 

6 c' 1 d 2 c 
b---cd6d~t = ° (K8) 
The solution of Euclidean Einstein equations for this metric is : 

a(t) = H^cosiHt), (E.9) 

1 2A 



b 2 [2(m + l)+n 2 (n 2 -l)]' 

c{&) = cos(0), 
,2 _ 2 



ni& 2 
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The static scale factor d of the torus T m is not fixed by these equations. 

One may consider S 1+ni xT 1 ^ 2 xT m type instantons in 10- D (l+rii + l + 2+m = 10) 
spacetime and see how the action varies. Using ( E.9| ) one can find the action for such an 
instanton. The action for an instanton S 1+ni x T 1 S ri2 x T m in D-spacetime dimensions is 
given by: 

e 1 4A t/ / 

&E = ^ttt x -^ x Volume 

16ttG d D-2 

4A 2vr( 2 + n i)/ 2 (F- 1 ) 1+ni 2^ 2+n2 )/ 2 (&f 2 ,2irA 1+m 



{—f +m (E.10) 



i6ttG d (d-2) r(24pi) r(^) v m s 

where ^4 is the size of the torii in string units. Using (|E.9| ), ( |3.4| ), (|3.3;), for the case D = 10 
and ri2 = 2 one obtains: 

-, 1 /7T \(3+m)/2 

S ^-(io^) nM „ i(ni/2) , (-„ l( 2 + ni )) (E.11) 

F. Locating the Preferred Inflationary State 

Consider the KKLMMT scenario. The Kahler potential is 

fC = -3ln[-i(p-p)]-ln[-i(T-f)]-ln(-if fi A Cl J . (F.l) 

The super symmetric minimum is given by D Z W = D T W = D p W = (where DiW = 
d{W + difCW), which lead to, respectively, 

£(!+]»,) -(K + tf^Jr.O (F.2) 

M£? - (ifz + KV)f = -,4e iap (F.3) 

M£? - (Kz + #V)r = ( -a(p - p) - 1 J ^e wp (F.4) 



Here we present solutions of Eqs( [F.2D , ( |F.3| ) and ( |F.4| ). To simplify the analysis we take 
r = ie~^ = i/g s , Pmin = i& while z and z' are taken to be real. Also we assume z'(z) ~ 0(1) 
and dz'(z)/dz ~ for small z. In the limit of physical interest where K/g s is large, we 
find from Eq.(jR2|) that 

In this limit, from the sum and difference of the remaining two equations, respectively, we 
find 



1 
9a 



aa = -3-LW(-l,-m) (F.6) 

- 3K'l'(o) A aa exp(-oo-) = f^j /?(m) (F.7) 

where we have introduced m = 3e _3 M|<7(0)A _1 | and /3(m) = LVKf-i -in) ^° ^ ave a 
solution, we need m < 3e -3 ~ 0.15. The function y(x) = — LW(— 1, — x) is the LambertW 
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functions, which is the solution of the equation ye y = x,x < e 1 , for y > 1. Using 
Eq.( [F.7| ) in z expression, Eq.(F.5), we find 



z = e -2n\g(O)\0(m)K/z'(O)K' _ (pg) 

We are interested in maximizing T = —Se — cMjVg in terms of fluxes, K,M and K'. 
Using the solutions found for z, p(a) and g s and using Eqs ( [4.12| ) and ( [4,14] ) for Se and 
cM^Vg , respectively, we find 

87r|6(0)|X/3(m) 4vr|e(0) \Kf3{m) 

T = Afi(m)e ^ cBf 2 (m)e k> (F.9) 



where 



B . *£ ( -£J)' \A\ia~f 

(2tt) 5 \3z'(0)J ' ' 



/i(m) = m 



3 (3 + LW(-l,-m)Y 

\LW(-l,-m)\ 3 



ff \- * \* + LW(-l,-m)ft , . 

t2(m) = m? 5 — (r-10) 

\LW(-l,-m)\2 



From the combination of -wwT = and -StjT = 0, we obtain 

3/i/ 2 - 2/i$ = (F.ll) 

pK/3(rn) 2 A him) 

e 2K> = — =±±± — '- 

3cBf 2 (m) 

The root of the first equation, independent of a and \A\ is m max ~ 0.06. It turns out 
that the condition -^rT = is not quite compatible with -ArJ" = 0. On the other hand, 
K' appears only in the denominator of the exponential terms in T ' ■ Since the fluxes are 
quantized we can choose the smallest possible integer value for K' to maximize J-\ that is, 
K' = l. 

From Eqs ( [F.ll ) and ( F.12 ), respectively, we find 



me 3 . A , ,„ 

M = -^r (F.12) 

3 '^(O) 1 V ; 

K _ 3 Uy/Zal him) e 3 A a A 

~ 47r0(O)/?(m) ^ (2vr) 4 / 2 (m) ' 3/(0) ' * cj " 

Correspondingly, the warp factor is given by 

[ R } [ (2vr)4 /sM'sz'W 1 6) { "- L6) 
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From Eqs ( @Alj ), ( |F.22[) and (|F.16| ), we find 

M P | 2 2e 6 m 2 /3(m^ 2 



M s 



i(2tt)V(0) 



3 + LW(-l,-m)|2a"2|A|" 



(F.14) 



To get a better understanding of the results, we will use the approximate values for (3(m), 
fi(m) and f2{m) at the maximum point m ma x ~ 0.06: 



P(m) ~ 0.3 



/i("i) 



10" 



/ 2 (m) ~ 5 x 10~ 3 



(F.15) 



Also the geometric quantities \G(0)\ and ^'(0) are basically calculable. We use the GKP 
approximation and set |^(0)| ~ -z'(O) ~ 1, so the expressions for M, K, the warp factor 
and Mpi/M s , respectively, simplify to 



M ~ 0.4 x \A\ 

K -0.8 x In (l0- 4 a 3 / 4 |A| 3 / 2 c" 1 

(^) - 115 x a" 3 / 8 | AP^c 1 / 2 
Mm 



(F.16) 



~3xlO- 4 a^ 3 / 4 |A| 



The value of function T at the maximum is 

t o sy m-is „-3/2 I 416 --2 

/mm ~ 2 x 11) a 'Ac 



(F.17) 



Finally, the COBE normalization (that is, the density perturbation 5h — 10 5 ) is given in 
Eq.(C.8) of KKLMMTas 



r a 



5h = C lN ^ (4^) 1/3 (?) 4/3 



M/ v i2' 



(F.18) 



where JV e is the number of e-folding and C\ ~ 0.4 in their model. Using Eqs (^8), ( 4.11| ) 
and ( |F.17 ) and taking N e ~ 60, we find 



5 H ~ 2 x W 7 a 1 ' 2 \A\~ 2 c 2 / 3 



(F.19) 



There is an interesting relation between T and 5h- One can easily show from Eqs ( [F.17] ) 
and ( |Fl9D that 



J~rr 



1.6 x 10 4 (57/ 



For 5h ~ 10 , we find 



^~ 10 



18 



(F.20) 
(F.21) 



This value for T is many order of magnitudes bigger than the value of J- for S 10 . 

To check the validity of the low energy supergravity limit, we must verify that g s < 1 
and r » l(in the unit of M s ). From Eqs ( |F.7| ) and ( F.16J ) we find 

1 sM 



g s 27rK'z'(0) 



0{m) ~0.1 x \A\ 



(F.22) 
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To satisfy the bound g s < 1, we take \A\ > 10. From Eq.( [F.6| ) we have 



cr 



— (3 + LW(-l,-m)) ~ — 

a a 



(F.23) 



where the approximation —3 — LW(— 1, — m) ~ 1.2 was used. We see that to have r >> 1, 
we need a « 1. For example the value chosen in KKLT for a = 0.1, results in r ~ 2. 

It is easy to generalize the above inflationary scenario to N pairs of D3 — D3 branes. 
We also add a very small Akklt like today's cosmological constant as is obtained in KKLT: 



A-> A = 87rGA/"T 3 



& 



+ Akklt • 



(F.24) 



For N 7^ 0, the first term clearly dominates and the second term may be ignored. This 
corresponds to A — > J\f~ 1 A and B — » N~ 5 ' 2 B in Eq.( |F.S| ). The extremum can be obtained 
exactly as before. M is given as in Eq.( [F.12" ) while K is given by 



K ^ K + 



InM . 



(F.25) 



8vrg(0)/3(m) 

So we see that the flux K is not very sensitive to the value of M. 

To compare the inflationary states with the KKLT vacua, we calculate J- for J\f = 0: 



Fkklt 



3tt 



2Gn Akklt 



1 



(2tt) 6 t 2 M 8 
2\/3 \/ Akklt. 



(F.26) 



To find an estimate for Tkklt, let take g s ~ 10 -1 , c ~ 10" 10 and A K kltIM 2 P[ ~ lO" 120 
corresponding to today's cosmological constant. One can easily verify that the T> term is 
exponentially larger than —Se, and 



Fkklt 



K M P[ 



(F.27) 



This means that the tunneling to a small A universe is severely suppressed, in contrast to 
that suggested by the Hartle-Hawking scenario. 
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